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Abstract 

We study the renormahzation of gauge invariant operators in large A^c QCD. 
We compute the complete matrix of anomalous dimensions to leading order in 
the 't Hooft coupling and study its eigenvalues. Thinking of the mixing matrix 
as the Hamiltonian of a generalized spin chain we find a large integrable sector 
consisting of purely gluonic operators constructed with self-dual field strengths 
and an arbitrary number of derivatives. This sector contains the true ground 
state of the spin chain and all the gapless excitations above it. The ground state 
is essentially the anti-ferromagnetic ground state of a XXXi spin chain and the 
excitations carry either a chiral spin quantum number with relativistic dispersion 
relation or an anti-chiral one with non-relativistic dispersion relation. 
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1 Introduction and summary 



The computation of anomalous dimensions of composite operators is of central impor- 
tance in many areas of physics, most notably condensed matter and particle physics. In 
particle physics, their phenomenological interest stems from hadronic processes at high 
energy. For instance, in the study of deep inelastic scattering P one is led to consider 
the operator product expansion of two hadronic currents. The anomalous dimensions of 
twist-two operators appearing in the expansion control the logarithmic deviations to the 
Bjorken scaling. Considering operators with low twist amounts to considering only oper- 
ators made up of very few fundamental fields F^y and q (two in most cases) and possibly 
many derivatives. These are the operators that are important in collider physics. 

On the other hand, the study of operators containing many fundamental fields has 
received new impetus in connection with the attempts to construct a string description 
of QCD at low energies W- There is probably still a long way to go in fulfilling this 
dream (if indeed the idea is viable), but for other gauge theories the answer has been 
found. We now have a string theory description of the maximally supersymmetric A/" = 4 
SU{Nc) theory in terms of type IIB strings in an AdSr, x background jHlll]- This was 
rather surprising at first since the J\f = 4 theory is not confining and there are no color 
flux tubes in the ordinary sense. The connection with string theory was instead obtained 
by identifying the string spectrum with the spectrum of anomalous dimensions of the 
theory jlj. 

After many checks performed in the supergravity regime, a new important step was 
taken in j3] where it was shown that by taking an appropriate limit of both sides of the 
correspondence one could match anomalous dimensions of near BPS operators with the 
corresponding string states. On the gauge theory side this amounts to considering long 
operators with large classical dimension and R-charge. These operators are constructed 
by taking the BPS operator ti Z'^ , where Z is one of the three chiral superfields making 
up the Af = 4 multiplet, and adding a few "impurities", i.e. replacing some of the Z 
with other fields in the theory. The analysis of near BPS states led to an agreement to 
all loops in the 't Hooft coupling^ asNc jH]- 

States with large quantum numbers should correspond to semiclassical string states 
whether they are near BPS or not. In fact, as noticed in [Zj one can get a qualita- 
tive agreement between the anomalous dimension of twist-two operators of large spin 
S computed at small agN^ and the energy of a classical rotating string in AdS space, 
valid for large values of asNc- In both cases, the anomalous dimension grows like logS", 
one of the most famous "trade- marks" of non-abelian gauge theories p. Of course, the 
dependence of the coefficient of log S on the 't Hooft coupling Os-ZVc is polynomial on the 
gauge theory side, as always in perturbation theory, whereas it grows like y/a^Nc in the 
semiclassical AdS analysis, a "screening" phenomenon that has already been observed 
e.g. in the analysis of the Wilson loops fH|. 

The above qualitative observations can be made quantitative by considering, once 
again, the A/" = 4 theory [HHTT]. In this case the string can instead be allowed to rotate 
in the compact dimensions. Giving the string two large angular momenta Ji and J2 

^We use the standard definition as = gyM/^""- 
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on the sphere (recall that the isometry group is 5*0(6) and has rank three) one obtains 
an expression for the energy that is analytic in Og^c and thus can be expanded around 
c^sNc = 0. In order to match this expansion with the perturbative analysis one needs to 
compute the anomalous dimension of operators that are far from the BPS condition. This 
is successfully accomplished by noticing that the matrix of planar anomalous dimensions 
can be interpreted as the Hamiltonian of an integrable spin chain - not only at one- 
loop [T^IIIH], but also at higher loops [H] - and its eigenvalues can thus be obtained by 
applying the Bethe ansatz [T^llinilinilini ■ For an introduction to this subject, c.f. [16j. 
(See also ^Tj for a computation of correlation functions from the Bethe ansatz.) 

Following this line of thought, Kruczenski showed that the agreement is already 
present at the level of the effective sigma model (see also [HI for related work). Namely, 
the effective action obtained by taking the long wavelength limit of the one-loop spin 
chain Hamiltonian and the action of the rotating string in the corresponding space-time 
background are the same. This leads to the exciting possibility of extracting information 
about the string dual from the perturbative study of anomalous dimensions and was the 
main motivation behind the previous work by some of us |20jj in the context of large Nc 
(non-supersymmetric) QCD. In [201 we showed that the closed spin chains constructed 
with only the selfdual field strength form an integrable subsector. In this paper we 
continue this investigations in the hope that it will reveal the nature of the excitations 
required in the description of the QCD string. 

Before moving to describing the results of our investigation, we must stress that the 
subject of integrability in the context of QCD has a long parallel history starting with 
the introduction and the study of operators on the light cone in ^J, (see also j23) and 
continuing with the discovery that their renormalization is described by an integrable, 
non-compact spin chain The early papers on the subject dealt only with "short" 

operators, i.e. operators containing a small number of elementary fields but with an 
arbitrary number of derivatives. We refer the reader to the reviews [26,27j for a survey of 
the literature on this subject. Very recently, the complete matrix of anomalous dimension 
for light-cone operators of arbitrary length has been given in [2H1 for theories with any 
number of supersymmetries. One of the results in the present paper is the complete one- 
loop matrix of anomalous dimension for large A^^c QCD without the restriction to the light 
cone. This was obtained by first solving the counting problem for operators and then 
truncating the known Af = 4 answer of jSHI but it can also be thought of as a "lifting" of 
the light cone results of j2Hj to the full conformal group. For earlier partial results that 
were useful in the initial stages of this investigation see jSO]. Recently, anisotropic spin 
chains have also made an appearance in this context |8T|l28]. 

The main results of the present work are the following: 

We give the complete one- loop matrix of anomalous dimension for large A^^c QCD 
and study its eigenvalues. We use a particularly convenient choice of basis for such 
operators that avoids the necessity of restricting oneself to light-cone operators. In 
order to use the intuition from condensed matter we will be thinking of the matrix of 
anomalous dimensions as the Hamiltonian of a spin chain. We are particularly interested 
in the "thermodynamic limit" consisting of studying operators with a large number of 
elementary fields. 

We are able to construct the exact (anti-ferromagnetic) ground state, corresponding 
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to the operator with the lowest (negative) anomalous dimension for any given length and 
we show that all gapless excitations above the ground state are described by an integrable 
system based on the conformal group SO{A, 2). There are in fact two equivalent ground 
states, corresponding to scalar operators constructed by only using the selfdual field 
strength or by only using the anti-self dual component. Operators constructed using 
both, or operators containing quarks, are separated by a "mass gap" from the ground 
state, i.e. their anomalous dimension remains higher than the ground state by a finite 
amount in the thermodynamic limit. 

With the above "stringy" motivation in mind we then analyze the structure of the 
operators that contain "gapless" excitations. These are the QCD equivalents of the BMN 
operators for A/" = 4. 

The analysis of the quantum numbers carried by the gapless excitations reveals some 
surprises. Focusing on the chiral ground state for definitiveness, we find, as already 
expected from [20], excitations carrying a (space-time) chiral index a corresponding to 
the "spinous" of the compact XXXi chain embedded in the conformal chain [SU{2)l € 
5*0(4,2)). As well known from the condensed matter literature, such excitations along 
the chain are characterized by a linear ( "relativistic" ) dispersion relation. 

But there is another type of excitation that was not accessible by the previous analysis 
which did not include derivatives. Namely, by "spreading" a covariant derivative along 
the ground state chain, very much like the impurities in the BMN case, we are able to 
construct new excitations along the chain that are characterized by a quadratic ("non- 
relativistic" ) dispersion relation. One would naively expect the quantum number carried 
by these excitations to be a space-time vector since a covariant derivative carries a space- 
time index /i. However, by writing /i = (a, a) and recalling that a propagates in an anti- 
ferromagnetic background, we are able to show that the presence of gives rise to two 
independent elementary excitations, one being the above mentioned spinon and the other 
carrying a space-time anti-chiral index a. It is well known since the work of Faddeev and 
Takhtajan that scattering off background spins can modify quantum numbers of 

excitations in anti-ferromagnetic spin chains. Here we found another manifestation of 
this phenomenon. 

The picture that emerges from the one-loop analysis has some intriguing similarities 
to the twistor string theory of Witten [HI] and it would be very interesting to make 
the connection more explicit. What we can say is that the thermodynamic limit of 
the one-loop QCD chain is described by chiral/anti-chiral twistorial excitations with 
linear /quadratic dispersion relation. Whether these are (part of) the correct degrees 
of freedom required to formulate the string description of QCD, what their space-time 
dynamics is and the connections (if any) with the work of [SI] remain a project for the 
future. 

The paper is organized as follows: 

We begin in Section |21 by discussing some basic preliminary facts about the renor- 
malization of composite operators, the simplifications occurring in the large A'^c limit and 
the definition of the spin chain Hamiltonian. 

In Section |S1 we discuss in detail the construction of the complete basis of operators 
that form the Hilbert space of our spin chain. Having such a complete basis makes the 
restriction to the light-cone no longer necessary. 
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Section ^ contains the complete one-loop chain of anomalous dimensions expressed 
in terms of projectors on irreducible representations of the conformal group. Recall that 
the use of conformal symmetry to classify the operator to one-loop is possible, even in a 
non-conformal field theory such as QCD, since the beta function is of order 0((as^c)^) 
whereas the anomalous dimensions are 0{asNc). The relevant mathematical formalism 
is summarized in Appendix El The coefficients of the projection operators are obtained 
by a truncation of the corresponding matrix in the Af = 4 theory with the additional 
modifications required by the different wave-function renormalization. We also give 
a table of anomalous dimensions of all conformal primary operators up to dimension 
seven for pure glue and five for those involving quarks. We conclude the section by 
showing briefiy how the same results could also be obtained by "lifting" the light-cone 
results of |2HI- To this purpose, we consider for definitiveness the gluon-quark coupling. 
Although this is not the way we arrived at the results, we feel it may be useful to include 
it in order to establish a connection with the previous literature. 

Section El presents the Bethe ansatz for the chiral sector of the theory. We show that 
the closed chain composed of selfdual gluon field strengths and derivatives is integrable 
and we give the Bethe equations corresponding to the conformal group S'0(4, 2). We 
begin our discussion of the possible excitations over the ground state. As a check of our 
computation we present all the Bethe roots corresponding to the conformal primaries 
in the integrable sector up to dimension seven. We conclude with some remarks on the 
open chain, which is also integrable for selfdual gluons in the bulk. We present the Bethe 
ansatz for the case of open chains without derivatives. 

Section IHl contains a more detailed analysis of the excitations over the ground state. 
First, we show that operators containing a mixture of selfdual and anti-self dual fields or 
quarks have an anomalous dimension that remains above the previously found ground 
state by a finite amount (gap) in the thermodynamic limit. This establishes that we have 
indeed found the true ground state. We then move on to study the gapless excitations. 
As mentioned before they are characterized by a dispersion relation that is linear or 
quadratic for chiral and anti-chiral objects respectively. 

In Appendix ^ we review the oscillator representation of the conformal group and 
the decompositions required in the paper. 

In Appendix IHl we make the comparison with the length two primaries some of which 
are well know from the literature on deep inelastic scattering. We do this partly as a 
check and partly as an illustration of how to use the Hamiltonian obtained in Section HI 

In Appendix we give a proof of integrability for the chiral sector of pure YM by 
solving the corresponding Yang-Baxter equations. 

2 Preliminaries 

In this paper we consider anomalous dimensions of local, gauge invariant operators for 
massless QCD: 

L = -l-tiF^ + iqpq, (2.1) 
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where Ff^i, = df^A^, — d^A^ — igYulA^, Ai\ is the gauge field strength, thought of as an 
iVc X iVc hermitian matrix, and q is the quark field, transforming in the fundamental 
representation of the gauge group SU{Nc) and possibly carrying a finite number of 
fiavors. It will be very convenient in the following to express such field content in the 
chiral basis by decomposing F^i, into its selfdual and anti-selfdual parts / and / 

F,. = a^JU + (2.2) 
and the quark field in its chiral and anti-chiral parts ip and x'- 

We employ the 't Hooft limit of letting the number of colors A^^c go to infinity and 
the usual QCD coupling qym to zero while holding the 't Hooft coupling asA^c finite. In 
fact, we will deal with one-loop perturbation theory in agN^ which must therefore be 
assumed to be small. 

The latter condition would make the limit rather simple were it not for the fact that 
we apply it to the study of anomalous dimensions of gauge invariant composite operators 
made out of a large number of elementary fields where even ordinary perturbation theory 
becomes quite involved.^ 

Let us briefiy recall how anomalous dimensions of composite operators are defined in 
quantum field theory. We use the conventions of throughout the paper. In particular, 
let us denote by (9 a multiplicatively renormalized operator, i.e. an operator that can 
be renormalized by multiplying it with a (divergent) renormalization constant. We set 

Cbare = ZoOren., (2.4) 

where Zo is constructed so that O^en. has finite correlation functions with the renormal- 
ized quark and gluon fields. The anomalous dimension of O is defined through 

d 

7o = Ai^log^o- (2.5) 
To one-loop in asN^ it takes the form 

7o = ^c?^^' (2.6) 

where co is a numerical constant to be determined. 

In general, operators carrying the same quantum numbers and having the same clas- 
sical dimension will mix under renormalization and we write: 

^bare = -^O^^ren.- (2-7) 



^Actually, with some more work it would be possible to obtain the mixing matrix at finite Nc but we 
decided not to include it in this paper as our main motivation is to learn about the stringy description 
where 1 /N^ plays the role of the string coupling. 
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The anomalous dimension also becomes a matrix and we write, in matrix notation: 

r = log Zo = Za'fi^Zo = ^H. (2.8) 

Multiplicatively renormalizable operators correspond to eigenvectors of F and their 
eigenvalues correspond to the anomalous dimensions: 

TO = -f^O. (2.9) 

The problem of the computation of anomalous dimensions is thus split into two parts: 
First one must derive the mixing matrix F and second one must find its eigenvalues. In 
the past few years new techniques, inspired by the study of TV = 4 SYM have been 
developed that bypass the computation of the divergences in the loop diagrams and 
allow a direct determination of F. These techniques are based on the analysis on the 
dilatation operator, thought of as one of the generators of the conformal group. The 
second part, involving the diagonalization of F, also requires analytical work since we 
are mostly interested in letting the size of the matrix go to infinity by considering "long" 
operators. Fortunately, even this second problem can be handled by analytic methods 
such as the Bethe ansatz for the cases of interest. 

We shall refer to a generic composite operator as a "chain" , where the name originates 
by the interpretation of such operators as a spin chain, as we will discuss in detail later. 

Neglecting for the moment the details of the Lorentz structure of these operators, 
such chains can be grossly divided into two main groups: (see Fig. H]) 




D'^q jjnp DPq 



Figure 1: Closed and open "spin chains" corresponding to purely gluonic operators or mesonic 
operators respectively. 

• Closed chains, i.e. operators schematically of the type 

tr(L)'^iF)...(D"^F), (2.10) 

• Open chains, i.e. operators schematically of the type 

(D"^g)(D"2F) . . . (D"^-iF)(D"^g), (2.11) 
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Notice that we will distinguish between different conformal descendants, i.e. will not 
identify operators differing by a total derivative. Therefore it does not suffice to ignore 
derivatives on the first site q as is commonly done in QCD. 

In the large limit the following three crucial simplifications occur: 

• If we consider two eigenvectors Oi and 02-, linear combinations of either type of 
operator described above, with anomalous dimensions and respectively, 

then V{0i02) = {lo^ +lo^){Oi02)- Hence, it is enough to consider operators that 
cannot be split into products of gauge singlets, commonly referred to as "single- 
trace" operators, although, in the presence of quarks the name "irreducible" would 
be more appropriate. 

• Open and closed chains do not mix. This can be seen by noticing that with the 
usual normalization, the mixing matrix between open and closed chains becomes 
upper triangular. Alternatively, one can show that by scaling the quark fields by 
an extra factor ^/Wl the mixing becomes block diagonal. We shall refer to the 
mixing matrix for open and closed operators as Fopen and Fdosed respectively. 

• Finally, the relevant Feynman diagrams for the computation of the anomalous di- 
mension are only those connecting neighboring fields (see Fig. |2)). Thus we can 
identify either type of basic operator with a chain of "spins", each spin corre- 
sponding to a basic building block D^f, D""ip, etc. The length L of the chain is 
conserved to one-loop and the matrix of anomalous dimensions can be regarded as 
a Hamiltonian of the spin chain with only nearest neighbor interactions. 




Figure 2: Renormalization of a four gluon composite operator (denoted by a cross) in double 
line notation. To the left is depicted a diagram connecting two adjacent gluon legs containing 
a closed color loop. To the right is shown a diagram connecting two non- adjacent gluon legs, 
not containing any color loop and thus subleading in 1/Nc- 

Thus we write, for the closed chain, with the usual periodic identification L + 1 1: 
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Tciosed - ^ (2-12) 
1=1 

and, for the open chain (L > 2): 



Topen = ( ^1^2 + X] + H^li^i ) . (2.13) 



1=2 



For completeness we should also add the rather trivial Hamiltonian for an open chain of 
length two 

Obviously, only ()2.12|) and ()2.13|) admit a continuum limit L oo. 

The matrix of anomalous dimensions is thus completely specified if we give the ex- 
pression for the various "link Hamiltonians" connecting two neighboring sites. 



3 Lorentz structure 

Before describing the link Hamiltonian, we must carefully define the space of operators 
on which it acts. This space must be constructed in such a way as to avoid double 
counting or missing some allowed operator. 

Seemingly different operators can be mapped into each other by repeatedly using 
either of the following two identities: 

[D^, D^] = -igYuF^.^ and D^^F^p + cychc = 0. (3.1) 

Moreover, we must identify operators that differ by the classical equations of motion: 

D^F;^ = -gYMqi.T''q and ^^D^q = Q. (3.2) 

Equations 1)3.11) and ()3.2j) will be used to maximize the number of building blocks in a 
given operator as much as possible in order to come to a unique basis. Before showing 
that, we should mention the reason why the equations of motions also need to be imposed. 
In [36j it is shown that if one separates gauge invariant operators into those not vanishing 
by the classical equations of motion, henceforth referred to as OS (on shell) operators, 
and those vanishing by the equations of motion (referred to as EOM operators), the full 
matrix of anomalous dimensions has the form:^ 

/pOS -pEOM\ 

^ ^ \ n pEOM • (3-3) 

V EOM/ 

Thus, although it might be necessary to add counterterms vanishing on shell in the 
renormalization of ordinary OS operators, we are assured that the former operators do 

■^To be precise, the EOM part might also contain BRST variations of gauge variant operators. 
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not alter the anomalous dimension of the latter since their mixing is upper triangular. 
From now on we shall consider only FqI. 

We now claim that a generic gauge singlet operator can be constructed, up to the 
above identifications, out of the following elementary building blocks: 

D"/ e V(„/2+i,„/2), and e V(„/2+i/2,n/2) (3.4) 

together with their complex conjugates. The fields /, ip and x are the chiral components 
of the gluon field strength and the quark fields discussed in the introduction. The 
labels {Si, S2) in ()3.4|) refer to the Lorentz quantum numbers and characterize a unique 
irreducible representation of the Lorentz group. The module 'V(s-i,S2) is a vector space on 
which the representation (5*1, 5*2) acts. In fact, these irreducible representations of the 
Lorentz group can be combined to form infinite-dimensional irreducible representations 
of the conformal group 50(2,4) by putting together all terms of the type, e.g. -D"/, 
n = ... 00, c.f. Appendix |X| The use of the conformal group is what will allow us to 
obtain the full matrix of anomalous dimensions to one-loop. 

To be specific, -D"/ written in the ordinary Lorentz basis would look like: 

= D^^, . . . D^^F^^^ - traces, (3.5) 

where F^^ is the selfdual component of the field strength. By "—traces" we mean that 
the tensor has been reduced to be totally traceless and the brackets around the indices 
correspond to total symmetrization. In the chiral basis, setting = cr'^'^Daa and 
= (J'^^fai3, one would write 

D"-/ = D^^ai ■ ■ ■ Da^a^fp-i + symmetrized, (3.6) 

where "-|-symmetrized" means that the tensor has been totally symmetrized in the un- 
dotted and dotted indices, respectively. One usually writes simply -D"' /«!... o„+2,ai... an 
for the resulting irreducible component in direct analogy with the representation labels 
(n/2-|-l,n/2), where each symmetrized spinor index contributes spin 1/2. 

To show the validity of the claim, let us look at a few examples. 

Consider first an operator containing the element D^F^p. We can decompose such an 
object into a totally anti-symmetric tensor, D[^Fj,p] vanishing by the Bianchi identity in 
fl3.1|) . a vector rj^yD^Fxp + rj^^pD^F^x proportional to the gauge current by the equations 
of motion, and thus quadratic in its field content, and a remaining tensor belonging to 
V(3/2,i/2) © V(i/2,3/2) and thus precisely of the type stated above. Hence: 

DpF,p = Df + Df + 0{qq). (3.7) 

Consider now an operator containing an element with two derivatives: DpD^Fpx. It 
is clear that the covariant derivatives can always be symmetrized, up to terms containing 
an extra F^^, due to ()3.1|) . Once the symmetrization has been done, the two indices of 
Fpx cannot further be anti-symmetrized with the ones of Dp or D^, because of fj3.1|) ^ and 
in total three indices will have to be symmetric. Finally, tracing any two indices will, by 

^This involves splitting off further F^^^s to bring Z?^ close to Fp\. 
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the equations of motion, yield a derivative of the current and thus a term quadratic in 
the fields. Perhaps the only slightly non-trivial identity occurs when the trace is between 
indices both belonging to the covariant derivative. In this case one must use once again 
the Bianchi identities to recover the equations of motion: 

= ^ + 0{qq) + 0{F^). (3.8) 

Again, the remaining tensors in the decomposition of the original tensor containing two 
covariant derivatives belong to the irreps (2, 1) © (1, 2). 

Similarly, the element D^q appearing at the end of an open chain will be equivalent 
to operators in V(i,i/2) © V(i/2,i) up to equations of motions. 

One can proceed by induction and show that: 

D^,... D^^F^p = D-f + D^f + Oiqq) + ^(F^), 

D,,... D„^q = D-^lj + D^x + 0{qqq) + 0{Fq). (3.9) 

Since ()3.9j] may raise but not decrease the length of a chain, one can proceed systemati- 
cally from the chains of lowest length (two) and show that the building blocks ()3.4j) are 
sufficient to construct uniquely all the on-shell operators. 

Note that we do not identify operators that differ by a total derivative. Although in 
some other applications this may be desirable, in our case we should keep them and we 
will identify them as descendants in the representation of the conformal group. 

4 Complete QCD spin chain 

We are now ready to discuss the form of the link Hamiltonian. We can think of it as a map 
from a generic pair of elementary building blocks previously discussed into two others. 
We can make use of Lorentz symmetry to simplify the form somewhat. The crucial 
insight is that the Hamiltonian can only map components from one irreducible module 
to the same component of a multiplet of the same type (or even the same multiplet). 
In general this is still a complicated problem because of the large number of Lorentz 
multiplets appearing in the product. Fortunately, for one-loop anomalous dimensions, 
conformal symmetry comes to the rescue by vastly reducing the number of occurring 
irreps as we shall now proceed to discuss. 

Although QCD is not a conformal theory, the use of the conformal group to classify 
composite operators in QCD has a long history (see [211 for a recent review). The 
conformal symmetry is especially useful at one loop because the /^-function, which is 
responsible for the breaking of conformal invariance, has leading term 0{{asNcY): 

Pia^N,) = f^§^^sN, = -y (4-1) 

and thus cannot affect the one- loop anomalous dimensions which are of order asNc- We 
shall make use of the conformal symmetry to write the full matrix of one-loop anoma- 
lous dimensions for large QCD. (The conformal algebra and its representations are 
reviewed in Appendix IX]). 
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Each fundamental field, namely /, ip and x and their conjugates, is a primary state^ 
of an irreducible module (or multiplet; a vector space transforming in some irreducible 
representation) of the conformal algebra. The remaining elements in the module (de- 
scendants) are obtained by acting with derivatives on the primaries. Thus, we have the 
following structure for the modules: 

{/01O2 ? ? f 01020^040102 • • • ) ; 

(Xoi) DXaia2ai^ D X01O203O1O2 • • • ) ) ('^•2) 

together with their complex conjugate modules. 

The product of two of the multiplets ()4.2j) can be decomposed into an infinite series 
of irreducible modules labeled by an extra integer j, the conformal spin. This is eas- 
ily understood in the oscillator representation reviewed in Appendix El We need the 
following decompositions (and their conjugates) 

j=-2 j=+2 



00 00 



j=-i j=+i 
00 00 

j=+l j=-l 

00 00 

Y'^ = Yj"^, ® = Yj^. (4.3) 

j=+i j=-i 

The explicit form of the irreducible modules appearing in the above decompositions 
is discussed in Appendix ^ The importance of these decompositions is that the one- 
loop link Hamiltonian commutes with the generators of the conformal algebra and thus 
assumes constant values on each of the irreducible modules. Moreover, each irreducible 
module appears only once within all the above decompositions and there will be no 
mixing problem between equivalent modules.^ 

Introducing the projectors Pj-^ etc. projecting on the respective modules we can write 
for each glue-glue link (not written explicitly) 

00 00 
H^F = Y (Pf + Pf ) + Yl + '''X)(Pf + Pf ), (4.4) 

j=-2 j=+2 

where the coefficients Ej need to be determined and I and X are the identity operator 



'^AU components of the Lorentz multiplets D'^ f, D^ip, . . . are conformal primaries. The conformal 
highest-weight state is distinguished by being also the highest-weight of the Lorentz multiplet. 

^This amounts to saying that for each irreducible module there is only a single (or 1x1 matrix of) 
coefficients instead of an n x n matrix in the case of an n-fold occurrence. 
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and the exchange operator, respectively: 



X 



I 



(D"/)(Z}™/) ^ {D^f){D-f). 



(4.5) 



Similarly, we write, for the quark-gluon and the gluon-quark link: 



oo 



oo 



H^F = ^ (Pf + Pf") + Ef\Pf + Pf ) 



j=~l j=+l 



oo 



oo 



= j2 (Pf + pf ) + Yl ^f^^? + pf ) 



(4.6) 



j=-i j=+i 



and for the length two quark-quark Hamiltonian: 



oo 



oo 



Y Ef{Pf + pf ) + J2 (Pf + pf ) 

j=-l j=+l 



(4.7) 



In writing the above equations we have already made use of various symmetries that 
allow us to identify the coefficients of the projectors on different modules. Perhaps one 
not so obvious symmetry is a "chiral" type of symmetry rotating / and / by opposite 
phases (see for a discussion in the context of A/" = 4 SYM). This symmetry is 
responsible for the closure of the chiral sector. We denote by A the quantum numbers 
of the fields under such symmetry transformation. 

The coefficients in equations ()4.4p4.fip4.7|l can be fixed in two ways, either by trun- 
cating the TV = 4 results of ^22j or by lifting the light cone results of P^25j . Many 
terms can also be explicitly computed or tested by considering chains of length two as 
discussed in Appendix |B] 

Before giving a derivation of the results we collect below the complete set of coeffi- 
cients of the two-site Hamiltonian, yielding the full matrix of anomalous dimensions to 
one-loop {h{j) = XlLi 1/^): 



£;f =2/i(j + 2)-f, 
Ep = Mj-2) + Mj + 2)-^, 

Ef = (-l)^' {h{j - 2) - Ah{j - 1) + Qh{j) - Ah{j + 1) + h{j + 2)) 
_ 6(-l)^'+^ 

" (j-i)j(i + i)(j + 2) 

Ef = h{j + 2) + h{j + l)-l, 
Ef = h{j + 2) + h{j-l)-l, 
Ef = 2hU + !)-§, 

Ef = h{j + 1) + h{j - 1) - |. (4.8) 



The spectrum of operators with the lowest classical dimensions can be found in Ta- 
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D L A 




[p, r, q] 




4 2 2 
4 2 2 
4 2 


(0,0) 
(2,0) 
(1,1) 


[0,- 4,0] 
[4,- 6,0] 
[2,- 6,2] 


11 + 

3 

7 + 

3 

0+ 


5 2 


(3/2, 3/2) 


[3,- 8,3] 


3+ 


6 3 3 
6 3 3 
6 3 3 
6 2 2 

6 3 1 
6 3 1 
6 3 1 
6 2 


(0,0) 
(1,0) 
(3,0) 

(3.1) 
(1,1) 
(0.1) 
(2,1) 
(2,2) 


[0,- 6,0] 
[2,- 7,0] 
[6,- 9,0] 
[6,-10,2] 
[2,- 8,2] 
[0,- 7,2] 
[4,- 9,2] 
[4,-10,4] 


1 + 

2 

3- 

2 
7- 
2 

14 + 
3 

7 + 

6 

11- 
6 

7- 

6 

21 + 
5 


7 3 3 
7 3 3 
7 3 1 
7 3 1 
7 3 1 
7 3 1 
7 2 


(3/2, 1/2) 
(5/2, 1/2) 
(5/2, 1/2) 
(1/2, 3/2) 
(3/2, 3/2) 
(5/2, 3/2) 
(5/2, 5/2) 


[3, - 9, 1] 
[5, -10, 1] 
[5, -10, 1] 
[1,- 9,3] 
[3, -10, 3] 
[5,-11,3] 
[5,-12,5] 


2+ and 2" 

17 + 
6 

17 + 

6 ^ 

12 ~12 

11+ and 1^" 
30 6 

197"^ and 
60 ^^^^ 12 
26 + 
5 



Table 1: A complete list of primary states for the purely gluonic sector of QCD up to = 7. The 

Dynkin labels \p, r, q] refer to the classical representation of the conformal group as explained 
in the text. The spin labels (5*1, 5*2) = {p/2, q/2) are redundant but we include them for clarity. 
D is the classical dimension, L the length of the operator and the chirality A counts the number 
of / minus the number of /. E is the anomalous dimension in units of . C denotes charge 
conjugation: C : F^y —F'^y In the language of spin chains it corresponds to reversing the 
orientation of the chain. For each state in the table with A ^ there is a corresponding state 
with opposite chirality, p and q exchanged and with the same anomalous dimension. 



D 


L 


A 




[P, r, q] 


E 


3 


2 


1 


(0,0) 


[0, -3,0] 


3 
2 


3 


2 


1 


(1,0) 


[2,-4,0] 


1 
2 


4 


2 


1 


(3/2, 1/2) 


[3, -6, 1] 


3 
2 


5 


2 


1 


(2,1) 


[4, -8, 2] 


13 
6 


5 


3 


2 


(0,0) 


[0,-5,0] 


4 

3 


5 


3 


2 


(1,0) 


[2,-6,0] 


-i and 1 


5 


3 


2 


(2,0) 


[4,-7,0] 


5 
3 


5 


3 





(0,1) 


[0,-6,2] 


1 
3 


5 


3 





(1,1) 


[2,-7,2] 


1 
3 



Table 2: A complete list of open chains with boundaries x • • • up to D = 5. 
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D L A 


(5*1, 5*2) 


r 1 

[p, r, q\ 


E 


3 2 


(1/2,1/2) 


[1,-4,1] 





4 2 


(1,1) 


[2,-6,2] 


4 
3 


5 2 
5 3 1 
5 3 1 


(3/2,3/2) 
(1/2,1/2) 
(3/2, 1/2) 


[3,-8,3] 
[1,-6,1] 
[3, -7, 1] 


25 
12 

1 
2 

1 



Table 3: A complete list of open chains with boundaries -0 ... -0 up to D = 5. For each state 
with A ^ there is a corresponding state with opposite chirality, p and q exchanged and with 
the same anomalous dimension. 

4.1 Reduction from SYM 

All the coefficients E^ ^ can be obtained from A/" = 4 SYM. The crucial observation is that 
the set of Feynman diagrams in A/" = 4 SYM encompasses all the Feynman diagrams of 
QCD. What is more, at the one-loop level, the additional propagating degrees of freedom 
of the supersymmetric theory appear only in a very restricted sense in the diagrams that 
are relevant to QCD. 

Let us explain this in more detail, starting with purely gluonic processes, i.e. inter- 
actions that couple only to gluons within a local operator and only emit gluons. The 
diagrams that contribute to the one-loop scaling dimension are of three basic types (see 
Fig. ED. 




Figure 3: Typical 't Hooft diagrams contributing to the one- loop scaling dimension of a four 
gluon operator. There are also "degenerate" diagrams where one of the vertices is connected 
directly to the gluon field appearing in a covariant derivative or a commutator inside the 
operator. 

Consider ffist the loop-interaction connecting to one field of the operator (Fig. El 
middle). In QCD, the particles in the loop can only be gluons (or the associated ghosts). 
Fundamental quarks are suppressed in the large Nc limit. In TV = 4 SYM, the gluons 
give precisely the same contribution, but also the adjoint scalars and fermions do couple. 
Hence we cannot simply read off the value of the gauge loop alone from the final result 
in A/" = 4 SYM. We merely know that it gives the same contribution to the scaling 
dimension for every gluon in the spin chain, but we shall leave the coefficient, E^ , to be 
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determined later. Of course it can be computed (and is well-known); at the end of this 
section we will show a way to derive it from conformal and flavor symmetry. To determine 
the contribution from the other two types of interactions (Fig. 01 left, right), we note 
that scalars and fermions in the vertices of A/" = 4 SYM always come in pairs. This 
means that for purely gluonic processes they can only contribute within internal loops 
of the diagram. These two diagrams do not have such loops, hence their contribution is 
exactly the same for A/" = 4 SYM as for QCD. 

We know that the complete one-loop dilatation generator of A/" = 4 SYM is given by 

oo 

H^=4 = ^^A^=4pA^=4 ^ith Ef=^ = 2h{j). (4.9) 

j=0 

The tricky part is the decomposition of the A/" = 4 invariant projector into projectors 
for the Af = conformal group. One finds that the module V^^^ contains only chiral 

combinations of two gluons Vj/, Vj/ with definite conformal spin j' = j — 2. We write 

pf =4 = . . . + p//^ + p//^ + . . . . (4.10) 

Consequently, the coefficient i^J-^ equals 

Ej^ = Ef=^^ + 2Ef = 2h{j + 2) + 2Ef (4.11) 

where E^ is the contribution from the missing loop-diagrams for scalars and fermions. 
For the non-chiral combination, the determination of Ej^ is more involved. The module 

V^^'* now contains the non-chiral combination Vj/ with conformal spin ranging from 
j' = j — 2toj' = j + 2. A second complication is that the two interacting particles may 
or may not change place. Finally, in A/" = 4 SYM they can transform into particles which 
are not part of QCD but carry the same quantum numbers of the conformal group. While 
these contributions are clearly dropped for QCD, they are essential for the definition of 
P^"^ as a projector. We would have to treat all two-particle states of A/" = 4 with the 
same quantum numbers as // to find the proper decomposition of projectors. This is 
rather involved and we just present the final result for the contributions Ej and Ej 
in (gSl). 

The situation for fundamental quarks requires yet another insight, because TV = 4 
SYM contains only adjoint fields. The point is that the gauge group structures at the 
one-loop level are very restricted. Again, they can only be of the three general forms 
depicted in Fig. El (forgetting for the moment the double line structure), where a vertex 
represents the structure constants A'* or f"''"^, depending on the type of line it attaches to 
(fundamental or adjoint). The loop diagram can again not be determined from Af = 4 
SYM and it gives rise to the unknown coefficient E"^. The other two structures can be 
distinguished by their symmetry under interchange of the particles. This is important, 
because in A/" = 4 the second structure can be transformed into the first one by means 
of a Jacobi identity. This means that for every contribution in the Af = 4 Hamiltonian 
we can derive the corresponding structure. This allows to derive also the terms for the 
fundamental fermions. Note that the Yukawa-coupling to the scalars in A/" = 4 involves 
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two different ffavors of fermions. It can therefore be suppressed by considering only one 
flavor of fermions. 

Finally, the unknown constants E-^" and E"^ can be obtained by demanding conserva- 
tion of the stress energy tensor and flavor currents. We find E-^" = and E"^ = — | in 
agreement with the /3-function of large A^c QCD and results of the next subsection. 

The power of the truncation from Af = 4 however shows up for integrability and the 
Bethe ansatz in Sec. 15.41 Here we needed to do detailed computations to obtain the 
energy coefficients E^^ , but the Bethe ansatz for QCD follows from the Bethe ansatz 
for A/" = 4 SYM straightforwardly. 

4.2 The lift from the hght cone 

The same results fj4.8|) can be reached by "lifting" the light cone results j22H2Sl- The 
state of the art for this technique is presented in j2Hl where the full expression for the 
matrix of anomalous dimensions for the so-called "quasi-partonic" operators is given. 
This amounts to restricting to the coUinear subgroup SL{2, R) of the full conformal 
group 5*0 (4, 2). Once the problem of classification and enumeration of the operators has 
been solved, their answer can be lifted to the full conformal group. 

Let us show how this works by lifting the expression for the gluon-quark terms which 
are the ones that cannot be fixed by looking at length-two gauge invariant operators. One 
works in light-cone coordinates -|-, — , _L where _L denotes the two transverse coordinates. 
The light cone decomposition is most easily obtained by introducing two light cone 
vectors n'^ and satisfying the Lorentz products n ■ n = n ■ n = Q and n ■ n = 1. 
One defines the -|- and — components by contracting with n and n respectively, e.g. 
= n^n'^F^y. (The notation is reviewed in [23123 •) 

The coUinear group is defined by the following three transformations, forming a sub- 
group of the full conformal group: 

X — 

x_ , x_ ^ x_ + c, x_ Ax_. (4.12) 

1 + 2ax- 

It is easy to show that this subgroup is isomorphic to SL{2,R). A generic field, either 
spinorial or tensorial, can be split into various components $ carrying a dimension d, a 
spin projection s (defined as the eigenvalue of the component S+_ of the spin operator) 
and a "coUinear conformal spin" defined as j = {d + s)/2, where we used the symbol j 
to distinguish the coUinear conformal spin from the conformal counterpart j previously 
used in the full S0{4, 2) context. We will see now that there is a close relation between 
these two quantum numbers. 

A certain field component $ can be evaluated on the light cone by introducing a real 
variable z: 

Hzn) = Y,^Dlm, (4.13) 

k=0 

where the Taylor expansion is a convenient way to keep track of the light-cone derivatives. 
A generic light-cone composite operator can thus be thought of as a polynomial in the 
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variables Zj, one for each elementary field. For instance, in [H^HITTj was considered the 
operator: 

S+{zi, Z2, Z3) = ^Y^q{zin){iF_L+{z2n) + Fj_+(z2ri)75)7+9(^3^), (4.14) 

for which the following Hamiltonian was given (changing slightly their notation to make 
the comparison with our formulas clearer): 

r = ^ [V(ii2) + f/F,(i23)] . (4.15) 

The coefficients Vqp and Upq, to be related to our E'^-l' and E^^ respectively, were found 
to be: 

V(j) = vl/(j + I) + vl/(j - I) - 2v[/(l) - I, 

Uf,0) = ^(i + i) + ^(i - I) - 2M/(1) - |. (4.16) 

In fl4.16|) we have denoted by \1/ the logarithmic derivative of the Gamma function "^{z) = 
r'{z)/r{z), related to the harmonic sum by: 

^(m) = /i(m - 1) - 7s, (4.17) 

'jE being the Euler constant that always cancels in the final expression. 

We now must relate the coUinear j to our previous conformal spin j. This can be done 
by observing that when acting on the rightmost link (i.e. for Upq) the primary fais'ip-y 
totally symmetrized has j = and contains the light cone component with j = 5/2 
[j = 3/2 from F and j = I from i/j). Thus in Upq we should set j = j ' + 5/2. Similarly, 
when acting on the leftmost link, (i.e. for Vqp) the relation should be j = j + 3/2 because 
now we have j = 1. 

Finally, in ()4.14|) there is an explicit factor of ^^ym in front of the operator. This 
means that the total anomalous dimensions are shifted by —11/6 in units of but, 
since the open chain has two links one must shift the contribution of each link ()4.16|1 by 
-11/12. 

Putting all this together yields back the previous results for E^-l" and E-^"^. In general, 
all the light-cone Hamiltonian can be lifted in a unique way to the full conformal group. 



5 Chiral operators and Bethe ansatz 

The conformal group provides a nice and compact bookkeeping tool for local operators 
in QCD which allowed us to write down the complete planar mixing matrix. The next 
obvious task is to diagonalize it. Though the complete diagonalization is beyond our 
reach, it will be possible to compute many interesting physical quantities. In particular, 
we compute the ground state of the mixing Hamiltonian to leading order^ in L and 
analyze the spectrum of small perturbations around it by the Bethe ansatz. Before 
going into the details of the general Bethe ansatz equations, we find it instructive to 
discuss a simple reduction of the full mixing problem. 

^In fact, one can systematically expand in powers of 1/L. The coefficients of and L^^ also have 
interesting physical meaning, being related to the boundary energy and the central charge respectively. 
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5.1 Chiral sector 



There are several types of operators that mix only among themselves and thus can be 
treated separately. An obvious example is a set of operators with the same quantum 
numbers such as spin, chirality or classical dimension. Though the length of the chain 
(the number of fields in an operator) is not a good quantum number, it is also conserved 
at the one-loop level, simply because producing an extra field requires an extra inter- 
action vertex and is thus suppressed by the QCD coupling. As in [201, shall first 
consider operators with minimal classical dimension for a given length. Such operators 
are composed of the gluon field strength and contain no derivatives: 

O = trF^^i,^ . . -F^^y^ , (5.1) 

A yet simpler subsector consists of chiral operators which contain only self-dual compo- 
nents of the field strength: 

= tTU,p,...f.^p,. (5.2) 

This relatively simple subset of operators contains the ground state of the full mixing 
Hamiltonian. 

Each entry faiPi in the operator represents one site of the chain of length L. Since we 
are excluding derivatives, the number of degrees of freedom per site is now finite. The 
three states faiPi form the spin-1 representation of SU{2)l- The spin operator is 

(SV)./^ = ^(crV^/ + ^/ay)A., (5.3) 

where cr*^'^, i = 1,2,3 are the ordinary Pauli matrices. The mixing matrix ()2.12j) acts 
pairwise on the adjacent sites of the chain. A pair of spins can be in the spin-0, spin-1 
or spin-2 state and the interaction Hamiltonian depends only on the total spin. From 
(I4.4l4.8jl we find that 

-11/6, if (S; + Sz+i)2 = 0, 
H,,,+i = { 1/6, if (S; + Sz+i)2 = 2, (5.4) 
7/6, if (S; + Sz+i)2 = 6. 

Using = 2, the Hamiltonian can be also written as [201 

Hz,/+i = |l + |S; • S^+i — |(Sz • Sz+i)^, (5.5) 

which is a Hamiltonian of the spin-1 quantum spin chain. Remarkably, this Hamiltonian 
is integrable ^4^ and its spectrum can be analyzed by the Bethe ansatz |HH111IS21ISS1 ■ 
It is convenient to use the basis in which the Sf^s are diagonal: 

/+ = /ll, /O = ^ (/l2 + /21) , /-=/22. (5.6) 



The operator 



On=tTf^ (5.7) 
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is an obvious eigenstate of the Hamiltonian with anomalous dimension 7^ = 

TagNcL /127r and spin S = L. This state is ferromagnetic (all spins aligned) and in 
the present case it is the state with the highest possible energy, which is easy to under- 
stand from ()5.5p : anti-alignment of nearest-neighbor spins lowers the energy, so the true 
ground state is anti-ferromagnetic. Quantum anti-ferromagnets are rather complicated 
systems even in one dimension, but in the present case the anti-ferromagnetic ground 
state can be found with the help of the Bethe ansatz PHIH^IHH] . 

It should be stressed at this point that integrability does not extend to the full 
Hamiltonian ()4.4|) due to the presence of the exchange operator X. Clearly, without 
such term, ()4.4|1 would be the sum of two independent integrable Hamiltonian, one for 
the chiral and one for the anti-chiral sector. The presence of X turns the system into a 
spin ladder by coupling the two sectors. However, it also spoils integrability as can be 
seen explicitly by constructing the candidate higher order conserved charges and showing 
that they do not commute with the Hamiltonian. In A/" = 4 SYM integrability is regained 
for the entire mixing matrix due to the presence of extra fields. In fact, we interpret the 
integrability of the chiral sector of QCD as a remnant of A/" = 4 integrability. 

5.2 The Bethe ansatz 

The Bethe ansatz describes all eigenstates of the spin chain in terms of elementary exci- 
tations around the ferromagnetic (pseudo)vacuum. The excitations close to the pseudo- 
vacuum (magnons) correspond to replacing some of the /+'s in ()5.7|1 by /o or /_. The 
magnons are thus created and annihilated by operators 

^ E a{l) ^ E ^"^'^'^"^ S+. (5.8) 

* m=l * m=l 

The momenta here are parameterized by rapidities: e*^*-'-* = (/ + i)/(/ — z), which is 
standard in the literature on Bethe ansatz. Each magnon reduces the spin by 1, so the 
state a^{li) . . . a"''(/A/)|0) has spin S = L — M. The rapidities k are in general complex 
because the magnons can form bound states with decaying wave function. The opera- 
tors ()5.8p create eigenstates of the Hamiltonian only asymptotically for very large chains 
and small number of magnons, when scattering of magnons on each other can be ne- 
glected. Nevertheless, the Hamiltonian can be diagonalized in a purely algebraic way if 
the creation and annihilation operators are appropriately deformed to take into account 
the scattering. The construction of the spectrum-generating operators, known as the 
algebraic Bethe ansatz, is rather involved and we refer to the original papers |521E] or 
to the review for the detailed derivation. Here we just quote the results. 

The rapidities of magnons should be all different and satisfy a set of algebraic equa- 
tions: 




which is basically the periodicity condition for the multi-magnon wave function. The 
right hand side of equation ()5.9p contains scattering phases. If scattering is neglected 
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the equation reduces to the quantization condition for the momenta: 

27rTi ■ 

p{lj) = 71 — 2arctan/j = — (5.10) 

Solutions of the Bethe equations {Ij, j = 1,...,M}, < M < L parameterize all 
eigenstates of the spin-chain Hamiltonian. The eigenvalues are 

There is a useful relationship between the momentum of a Bethe root and its contribution 
to the anomalous dimension: 

lil) = ^p'il)- (5.12) 

This relationship holds for arbitrary compounds of Bethe roots and will be very useful 
in calculations. 

The Bethe states which correspond to QCD operators should satisfy an extra condi- 
tion 

n^=i. (5.13) 

j=l 3 

which guarantees that the state has zero total momentum and hence is translation in- 
variant. This condition reflects the cyclicity of the trace in the operators ()5.2j) . 



5.3 Ant i- ferromagnetic ground state 

We are now in a position to construct the true, anti-ferromagnetic vacuum. Again we 
omit many details that can be found in the original papers |1^44j . The ground state has 
spin zero and thus contains L magnons. All their rapidities should be different, so the 
construction of the anti-ferromagnetic vacuum of the spin chain is similar to the filling of 
a Fermi sea. As we mentioned earlier, magnons with the same momentum form a bound 
state, so their rapidities become complex. The left hand side of the Bethe equation ()5.9j) 
then is not a pure phase and in the thermodynamic limit L — » oo it either blows up or 
goes to zero. This has to be compensated by a zero or a pole in the scattering amplitude 
on the right hand side. A pole or a zero arises when some of the Bethe roots differ 
by i. The bound state of k magnons is thus described in the thermodynamic limit by 
an array of Bethe roots with a common real part and integer or half-integer imaginary 
parts. Such a compound of Bethe roots is usually called a fc-string. The 2-strings play 
the most important role for the spin-1 chain. They are pairs of roots at 

hj-i = \j + V2, hj = Xj - V2, (5.14) 

where the centers of the strings Xj are real numbers. The ground state of the spin-1 
chain is the Fermi sea of 2-strings 
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To find the distribution of tlie 2-strings, we first multiply the Bethe equations for 
both roots and rewrite them in the logarithmic form: 

, , , 27™,; 1 v-^ , , , , X .X 

P2(A,) = + ;^ E'^2,2(A,- - Afc), (5.15) 



where P2(A) is the momentum of the 2-string: 

7 \ -\- .q,- /9 

(5.16) 



ir,,(x) A + z/2 A + 3^/2 
A - i/2 A - 3i/2 

which is equal to 



2A 

P2(A) = 27r — 2 arctan2A — 2 arctan — . (5-17) 

The branch of the arctangent is chosen such that — 7r/2 < arctan A < 7r/2. The phase 
ambiguity in choosing the branch is reflected in the mode numbers rij which must be 
different for different strings. For the scattering phase of the two strings we find after 
some calculations: ^ 

(^2, 2(A) = Svr — 4 arctan A — 2 arctan — . (5.18) 

Since in the vacuum all available one-particle states are occupied, we can put rij = j, 
j = l,...,L/2 in (j05|l . 

In the thermodynamic limit, the distribution of Bethe strings can be characterized 
by a continuous function A = X{x) of the variable x = j/L or by the density 



A'(x) 



(5.19) 



The thermodynamic limit of the Bethe equations for rapidities is an integral equation 
for the density: 



11 3 1 _ 

2 A^TTTi + 2 A^T974 " ""^^ ^ 



+00 



L(A-0' + l (A-0' + 4 



(5.20) 



This equation is derived by subtracting ()5.15|) for the (j + l)-th string from the equation 
for the j-th string and taking the difference Xj — Xj+i to zero. The equation can be easily 
solved by the Fourier transform: 

/•+°° dq e-'^^ _ 1 
P^^^ ~ 2^ 2^^ " 2cosh7rA ' ^^'^^^ 

The density is normalized as 

dXpi\) = l, ^5-22) 
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so that the ground state contains precisely L/2 2-strings or L roots and is therefore a 
spin zero state. The ground-state energy is, to leading order in L: 

which is the lowest possible anomalous dimension for operators of a given length. In 
the subsector without derivatives, the complete Hamiltonian is also bounded from above 
by the ferromagnetic vacuum. However, this is an artifact of the truncation - including 
derivatives it is possible to arbitrarily raise the anomalous dimension of operators of a 
given length. For instance, in the case of twist-two operators anomalous dimensions grow 
like log S where S" is the number of derivatives. 
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(5.23) 



5.4 Operators with derivatives 

We now consider adding covariant derivatives to the chiral operators discussed above. 
The resulting set of operators turns out to be the largest integrable sector in the pure 
Yang-Mills theory. At the same time this sector contains all low-energy excitations 
around the anti-ferromagnetic vacuum discussed in the previous section. We will show 
below that the other, non-chiral modes are separated from the vacuum by a gap. The 
most general operators in the chiral sector have the form 

O = tT{D"''f)...{D"'^f). (5.24) 

This sector is closed under renormalization at one-loop and integrability follows from 
integrability in A/" = 4 SYM by a simple argument: The operators ()5.24|) also form 
a one-loop closed subsector of = 4 SYM.* Therefore not only the mixing matrix for 
()5.24j) is inherited (up to a constant shift), but also its integrability.^ The mixing matrix 
can now be diagonalized by a Bethe ansatz. As usual, the Bethe equations are completely 
fixed by group theory and can be read off from the general result of |45'j, see also [1^ . 

The mixing of operators ()5.24|) is described by an 5*0(4,2) spin chain with spins 
in the representation whose Dynkin labels are \p,r,q] = [2,-3,0]. The left and right 
Dynkin indices refer to twice the Lorentz spins p = 2Si and q = 2S2, whereas the central 
Dynkin label is given by r = —D — Si — 82- The spectrum of the integrable spin chain 
with these symmetries is characterized by three sets of Bethe roots Uj, Ij and rj. The 

^Similar ideas have been pursued by R. Argurio (private communication). 
^In Appendix ICl we present an independent proof of this statement. 



22 



+2 —3 +2 —3 

Ml Mu Mr Ml Mu Mr 

Figure 4: Reduction of the distinguished Dynkin diagram of 5[/(4|2,2) to SU{2,2). Above 
the diagrams we have indicated the Dynkin labels of the spin representation and below the 
excitation numbers of each simple root. 

Bethe equations can be inferred from 



1 , ■\ L Ml J Mu J ■ /n 

_ -TT - /fc + ^ -TT ~ ""fc - V 2 

Uj — 3i/2\^ ^ Uj — Uk + i Yj Uj — Ik — i/2 -pr uj — — i/2 
Uj + 3i/2 J J-j^ Uj -Uk-i Uj - lk + i/2 Uj -rk + i/2' 

Mr , ■ Mu ■ /o 

nr~ — Tu + I -w-r — ui; — I 2 , , 
^^n^ ^- 5.25 
-n-t -Uk + 1/2 

These are precisely the Bethe equations of A/" = 4 SYM when one truncates the 
supergroup SU (2, 2|4) down to the bosonic part SU (2, 2). For this purpose one needs to 
consider the 'Beast' form in ^H] and remove (or simply not excite) the fermionic node 
along with the three nodes of the internal symmetry SU (4) from the distinguished Dynkin 
diagram of SU{2,2\4), see Fig. |31 The remaining Cartan matrix for the construction of 
the Bethe ansatz is the one of SU{2, 2) and the spin labels are [2, —3, 0], i.e. ()5.25|) . 

If one reduces further by not exciting the u and r roots, one returns to the Bethe 
equations ()5.25p for the SU{2)i sector. Equivalently, removing I and r roots leads to the 
SL{2) Bethe ansatz reviewed in PHIEZ!- 

The derivation of the momentum constraint and of the anomalous dimension follows 
exactly the same route as in the previous section. We simply present the results, which 
agree with A/" = 4 SYM up to an overall energy shift proportional to L. The cyclic 
states are in addition subject to the momentum constraint: 

The anomalous dimension is 




The Bethe state with Mi /-roots, w-roots and Mr r-roots is the highest weight in 
the 5*0(4,2) representation with Dynkin labels 

[2L + M„ - 2Mu -3L + Mi + Mr - 2M„, M„ - 2M^]. (5.28) 
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D L 


(5*1, 5*2) 


[p, r, g] 




roots 


4 2 


(0,0) 


[0,-4,0] 


-11/3 


/l,2 = ±z/V3 


4 2 


(2,0) 


[4, -6,0] 


7/3 


(no roots) 


6 3 


(0,0) 


[0,-6,0] 


1/2 


/i = 0, /2,3 = ±^ (singular state) 


6 3 


(1,0) 


[2,-7,0] 


-3/2 


/l,2 = ±VV5 


6 3 


(3,0) 


[6,-9,0] 


7/2 


(no roots) 


9 

U 


r3 ^^ 


[U, iU, ZJ 


14/3 


7/1 o I-3/a/2S 


7 3 


yoj z, ±/ 


f3 —9 11 

[O, J, ±J 


2 2 












/ _ 7 _ -r , / l+jV399 


7 3 


(5/2, 1/2) 


[5,-10,1] 


17/6 


Ml = /i = 



Table 4: Anomalous dimensions (7 = ^^:^E) and Bethe roots for the lowest dimensional chiral 
closed chains (^ = L in Table 

In other words it has 

D = 2L + Mu, Si = L + - Mu S2 = \Mu - M^. (5.29) 

The reference state with no excitations (the pseudo- vacuum) corresponds to the operator 
fj5.7|) that was identified in Section El as the ferromagnetic vacuum: 

On = tr//i=tr/^. (5.30) 

Adding a w-root corresponds to adding a derivative D^i to the above operator. An /- 
roots flips one left spin: 1 — 2, and an r-root flips one right spin: 1—^2. Since the 
Lorentz spins in ()5.29j) cannot be negative, the numbers of roots of different types are 
constrained by 

Mi<L + \M^, Mr<\M^. (5.31) 

Keeping only Z-roots brings us back to the Bethe equations of the XXXi chain discussed 
in [in. 

We can reproduce from the Bethe equations all anomalous dimensions of the operators 
in Tabled that belong to the chiral sector, i.e. those satisfying A = L. The results are 
given in Tabled for comparison. 

Adding the derivatives does not change the ground state, since w-roots give positive 
contribution to the energy. In Section IHl we will show that adding / to a chiral operator 
raises the energy by substantially larger amount. Thus, the ground state and all the 
low-energy modes of the mixing matrix are described by an integrable system. 

5.5 Open chains 

Before moving on to a more detailed discussion of the spectrum of the closed chain, we 
would like to make some comments on the open chain. (For earlier work on operators 
of this type we refer the reader to [231231111 •) The open chains with quarks at the ends 



24 



are also integrable if the gluon operators in the middle are chiral (see Appendix E])- Let 
us consider operators that have no derivatives: 

O = g/aaft • • • fai^-ih-iQ- (5-32) 
The mixing matrix for these operators is described by the following open spin-1 chain: 



^ -^ + Si-S2 + S,_i-s, + }^ 

l> 1=2 



7 1 1 



(5.33) 

where Lr = 0, 1, 2 counts the number of anti-chiral quarks. Here s depends on the 
chirality of the quark ()2.3|) that is inserted at each end of the spin chain: For chiral 
quarks we set s = cr/2, while anti-chiral ones have no SU{2)i^ spin and thus s = 0. 
Both types of boundary interaction are integrable. In the case with two anti-chiral 
quarks, we can take either the spin-0 or spin-1 combination of SU{2)ji without affecting 
the anomalous dimension. The chain with the boundary interaction has been studied 
recently in jlB] (see also [IH] for a discussion of open spin chains in a supersymmetric 
context). 

The Bethe ansatz for the system described in ()5.33|) is 

, , -N 2L-2~Lr M J 7 I • 7 I 7 I • 
k=l 



Ij I J 1^—1 ^ ~^ ^ 



and its eigenvalues are given by 



a,Nr I 7L-ALr - 11 ^' 



7 



2. \ 6 ^fTl)' f^'^^) 



k=l 



There is no momentum constraint. 

We can reproduce from the Bethe equations all anomalous dimensions of the operators 
in Table EEl that belong to the chiral sector, i.e. those satisfying A = L — 1 — L^. The 
results are given in Table |S1 for comparison. 

We will not study the thermodynamic limit of the open spin chain here, but we 
should mention that the anti-ferromagnetic phase of open spin chains exhibit rather 
unusual behavior under certain circumstances 



6 Spectrum of excitations 

We can now study the spectrum of excitations around the anti-ferromagnetic ground 
state of the spin chain. In part, this can be done using techniques that are well known 
from the literature on one-dimensional anti-ferromagnets. Let us first summarize our 
findings: There are three types of excitations^^ (Fig. EI): 

^"^In this context, an excitation refers to a distortion of the wave-function which carries an indepen- 
dent momentum. Physical states composed from such excitations are subject to constraints: A single 
excitation may be unphysical, see below. 
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D L Lr 




[P, q] 


E 


roots 


5 3 
5 3 
5 3 
5 3 


(0,0) 
(1,0) 
(1,0) 
(2,0) 


[0,-5,0] 
[2,-6,0] 
[2,-6,0] 
[4,-7,0] 


-4/3 
-1/3 
2/3 
5/3 


/i = 
h = l 
(no roots) 


5 3 1 
5 3 1 


(1/2,1/2) 
(3/2, 1/2) 


[1,-6,1] 
[3,-7,1] 


-1/2 
1 


h = 1/v^ 
(no roots) 


5 3 2 
5 3 2 


(1,0) 
(1,1) 


[2,-6,0] 
[2,-7,2] 


1/3 
1/3 


(no roots) 
(no roots) 



Table 5: Anomalous dimensions (7 = ^^^^E) and Bethe roots for the lowest dimensional open 
chiral chains {A = L — 1 — Lj. in Table 015]) . 

• Excitations with spin quantum numbers (0, 1) separated from the vacuum by a 
finite gap: e ~ const. These excitations correspond to changing a / into a /. 

• Excitations with spin quantum numbers (1/2, 0) and relativistic dispersion relation: 
e p 1/L. These modes are the low-energy states of the anti-ferromagnetic 
XXXi spin chain which are usually called spinous or spin waves. A single spinon has 
fractional spin and can therefore never appear on its own P^I42H44j . Furthermore, 
it obeys a non-standard exchange statistics. 

• Excitations with spin quantum numbers (0,1/2) and non- relativistic dispersion 
relation: e ~ ~ 1/L^. These modes have the lowest possible energy ~ 1/L^ 
and are turned on by the insertion of a covariant derivative into the operator. The 
derivative has both a left (chiral) and a right (anti-chiral) index but, since the left 
spin propagates in the anti-ferromagnetic vacuum, the left and right excitations 
propagate independently. The left excitation dissolves into the ground state and 
thus the only new element is the right excitation of classical scaling dimension 
1. These excitations are somewhat analogous to magnons in a ferromagnet and 
can form multi-particle bound states. They are probably the most surprising and 
potentially interesting result of our analysis. 

Before we move on to discuss the nature of each excitation in the next subsections, 
it must be stressed that, due to the highly non-trivial nature of the anti- ferromagnetic 
vacuum and to the zero momentum constraint, not all combinations of excitations are 
allowed. Only certain combinations of the elementary excitations above correspond to 
physical states that can be identified with QCD operators. This fact is already known in 
the case of spinous from the literature on anti-ferromagnets and is responsible for their 
unusual statistic We shall see more examples of this fact. 

6.1 Adding an anti-chiral impurity 

We start by showing that replacing one / by an / in the chiral spin chain gives rise to 
a mass gap, i.e. the chains containing such "impurities" have an energy that is larger 
than the vacuum by a finite amount in the thermodynamic limit L — > 00. For simplicity. 
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anomalous 
dimensions 



open chains and anti-chiral impurity ~ L° (gap) 



-==- spinous ~ L ^ 



gapless excitations 



IT derivatives ~ L ^ 

anti-ferromagnetic ground state 



Figure 5: Spectrum of anomalous dimensions constructed over the chiral anti-ferromagnetic 
ground state. An identical tower of excitations can be constructed over the anti-chiral anti- 
ferromagnetic ground state. 



we will work in the absence of derivatives, but it will be obvious from the discussion in 
Sec. 16.31 that the presence of derivatives can only increase the gap. 

For later use, recall that anti-ferromagnetic closed and open spin chains (with "free" 
boundary conditions) of length L have a ground state energy behaving as |52j : 

^ciosed(Iv) = aL + 0(1/L), Eopcn(i^) = aL + h + 0{l/L). (6.1) 

The constant a depends on the constant shift in the Hamiltonian. With our normalization 
it has been fixed in ()5.23|) to be a = —^^7^ < 0- '^^^ constant 6 > is the surface energy 
(not present for the closed chain with periodic boundary conditions). The terms of order 
0(1/L) are also of physical interest in many contexts, being related to the central charge 
of the system, but they will not be needed in the following arguments. 

In Sectional we have shown how to express the Hamiltonian in the chiral sector as a 
spin chain. In order to study the effect of impurities, we must now write the Hamiltonian 
containing both / and / (but no derivatives) as a spin chain (or better, a spin ladder). 
This is the Hamiltonian found in [2QJ to describe the renormalization of operators without 
derivatives. This Hamiltonian is, of course, the restriction of our present Hamiltonian 
()2.12j) to the subsector without derivatives. 

The spin system can be represented as a spin ladder. We have already described the 
action of S on / in ()5.3|) . The same operator acts on / in the obvious way. But now 
we also need to describe operators switching / with /. This can be done by introducing 
another independent spin operator r corresponding to a spin 1/2 representation and 
acting on each component of / and / (thought of as a doublet) as: 

tV = /, t+/ = 0, t-/ = /, 

r'f = -f, r+f = f, T-/ = 0. (6.2) 
The two spins ti and S/ can be visualized as sitting on a spin ladder, two parallel spin 
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chains connected by links at each site. The mixing matrix can be written as 

r = ^ E {5 (I + -hl^) + 5 (I - ^^n, (6.3) 

where h; ^+1 is the integrable spin one hnk Hamiltonian ()5.5|1 found in the previous 
section: 

= [fl + ^S, ■ S^+i - i(Si • Si+i)2] . (6.4) 

For chiral states, such as O in ()5.2|) . t'I^O = 0, rfO = O, and the mixing matrix reduces 
to 

Denote for convenience by {A = +, 0, — ) the spin one triplet introduced in 
and write the full wave function as a linear combination 

L 

^ = ^ ^M-A„.,BA„+^...A^ ^^^jA, _ _ _ jA„.,jBjA„+^ . . . /^^), (6.5) 



n=l 



in terms of some coefficients ipn- Clearly, one could always bring the impurity to the 
beginning of the trace, but it is convenient to have it at an arbitrary position. 
We want to solve the eigenvalue problem 

= Ei„ityiL)^, (6.6) 

where £'impurity(-f') denotes the energy of the chiral spin chain with L sites, one of which 
is replaced by an impurity /. Decompose 

tr(/^i . . . /^n-i/^/^n+i __jA,^^J2 tr(/^Hf,) (6.7) 

n 



where Q and are respectively the eigenvalues and eigenvectors of the open XXXi spin 
chain with L — 1 sites and (j^"+^---^LAi-An-i ^j^^ overlap coefficients. 
Partially projecting the original coefficients tpn onto the new basis: 

and doing some index manipulations we get the secular equation 

+ Ei^p,rity(^))xf = lixti + X^+i) (6.9) 

which has the usual plane wave solution. Of course, one must pick a direction in the 
anti-chiral index, say B = 0, and the others are obtained by SU{2)ji rotation. 
Thus, setting Xn = ^^^^ get the dispersion relation 

-Eimpurity(-^^) = f2 + 3 (1 - COS p) . (6.10) 

The energy of the open chain on the RHS must be evaluated for L — 1 sites because 
removing the impurity not only opens the chain but also lowers the number of sites by 
one. The lowest value for the energy is thus reached for p = 0. The important point 
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is that we can now use the bounds (c.f. (jfi.lll ) together with the fact that a < < 6 to 
show that 

^impurity(i^) = ^closcd(^) + (6.11) 

thus estabhshing the presence of a mass gap for these excitations. 

Let us move on to consider now the gapless excitations, all to be found in the inte- 
grable chiral sector. 



6.2 Spinons 



Here we review the spectrum of excitations around the anti-ferromagnetic vacuum of 
the XXXi model |12]. The excitations correspond to creating holes in the Fermi sea 
by removing one or more 2-strings^^ from the vacuum and/or to adding fc-strings with 
k ^ 2. They obey a non-trivial exchange statistics inHEHl- 

We have already obtained the density of 2-strings p(A) in the ground state in ()5.21|) . 
If one of the 2-strings is removed, then instead of ()5.2()|1 the density of 2-strings satisfies 



2 A2 + 1/4 
+2 



+ 



2 A2 + 9/4 

C?^Phole(0 

The density of 2-strings now has the form^^ 



7rphoic(A) + — 5(A - /i) 



+ 



(A-O^ + l (A-0' + 4, 



Phoio(A) = p(A) + - a-hoie(A 
and can again be found by the Fourier transform: 



P) 



(6.12) 



(6.13) 



0"holc(A) 



dq 



2ti 4cosh2f 



-5(A) + 



A 



2-i\ 



- m 



2 sinhyrA 

1 + zA 



- m 



1 -i\ 



(6.14) 



where ^I/(A) = r'(A)/r(A). The energy and the momentum of the hole can be readily 
calculated: 

2AA 



Phole(p) 



+00 



d\ o"hoic(A — p) ( vr — arctan2A — arctan ~^ j 



TT 



— arctan sinh tt/x. 



^holc(A') — — 



(iA (Thole (A - P) 



+ 



A2 + 1/4 A2 + 9/4y 2cosh7r/i 



TT 



(6.15) 



^^It is instructive not to think of holes as removed 2-strings, but rather as gaps in the sequence of 
mode numbers (e.g. nj = j, n^+i = j + 2) which was defined below IjS.f 8|) for the ground state. This 
explains why the hole excitation carries spin —1/2 and not —2. 

-'^^In order to avoid confusion, we stress that we will always be concerned with the density of 2-strings. 
By the notation phoic we mean the density of two strings in the presence of a hole in the sea of two 
strings. Similarly for other type of impurities. 

^•^The integrals are most easily done by using the Fourier representation for the density and the fact 

t^^at pUe(A*) = -2£holo(M)- 
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Thus the dispersion relation, 

echoic (p) = - sinp, (6.16) 

is hnear for the low-energy states. 

We can also compute the quantum numbers of the hole. And here we encounter a 
surprise |32j. Each individual 2-string contains two Bethe roots and thus has spin 2, but 
removing a 2-string from the Fermi sea distorts the distribution of the other strings and 
this distortion partially screens the spin of the hole! Naively, since we have removed a 
2-string, we would expect the correction to the density (Jhoio to be normalized to minus 
one. Instead we find: 

1 

c?A(Thoic(A) = -- , (6.17) 

so the perturbed density contains not L — 2, but L — 1/2 Bethe roots. Of course a state 
with a fractional number of roots or even a fractional number of 2-strings does not make 
sense and the holes in the sea of 2-strings can only be created in sets of four.^^ Yet each 
hole is an independent excitation of the anti-ferromagnetic spin chain, which has the 
dispersion relation (j6.16j) and spin 1/2 [22] • The Lorentz quantum numbers of the hole 
are (1/2,0). 

We now want to investigate what happens if, instead of making a hole, we add one k- 
string with k ^ 2. Of course, since the number of /-roots is already maximized in the anti- 
ferromagnetic vacuum, this by itself will not lead to a physical state but the modification 
in the density of 2-strings can be investigated independently of this issue to order 0{1/ L). 
If a /c-string with 7^ 2 is added to a distribution of the 2-string, the equation for the 
density of 2-strings takes the form (now with Pkstvi.^) = p(A) + \ o"fc_str(A — jj)): 

11 3 1 

+ o x9 , n/. = 7rpfe_str(A) 



2 A2 + 1/4 2 A2 + 9/4 
k-2 1 k 1 k+2 1 



2L (A-^)2 + (A;-2)V4 L(A-^)2 + A;2/4 2L (A - ^)2 + (A; + 2)2/4 



/+00 
■00 



+ 



[(A-O^ + l (A-0' + 4j ■ 



(6.18) 



The middle line comes from the scattering of a 2-string on the fc-string whose rapidity fi 
should be determined self-consistently and depends on the rapidities of other excitations. 
The middle line can be obtained by taking the derivative of the scattering phase for this 
process: 

2A 2A 2A 

^2kM = '271(2 ~ Ski) — 2 arctan 4 arctan 2 arctan . (6.19) 

k — 2 k k + 2 

For the density perturbation we find: 



^''We have tacitly assumed that L is even. For a spin chain with odd length L, the filled Fermi sea is 
not physical as it contains a half- integer number L /2 2-strings. Here, two holes would lead to a physical 
state. 
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type 


Ml Mu M2 


Si S2 D 


e{p) 


hole 

1-string of /'s 
A;-string of /'s, n > 2 


-1/2 1/4 
1/2 
k-2 1 


1/2 

-k + 2 


~7rb|/2 


1-string of m's 
fc-string of m's, n > 1 


1/2 1 -1/4 
1 k -1/2 


1/2 1 
A;/2-l k/2 k 


~ p^/2k 



Table 6: Chiral excitations and their quantum numbers. M2 indicates the number of the 
2-strings removed from the Fermi sea. 



and 

""^--'^^^^ 2n + (A; - 2)2/4 ' 

for k > 2. The modification in the density of 2-string in the presence of one such A;-string 
is thus normahzed as 

1 



rfAai_,tr(A) = --, (6.22) 
J -00 ^ 

and 

/"+00 

rfAafc_,tr(A) = -1, (6.23) 



so that effectively the 1-string adds no roots (—1/2 of a 2-string plus the single Bethe 
root) and the /c-string adds k — 2 roots to a Bethe state. However, all fc-strings do not 
contribute to the energy and momentum and thus do not correspond to propagating 
degrees of freedom. They are responsible for the spin and internal degrees of freedom of 
the magnons. 

More generally, a physical 2A^-magnon state will contain 2A^ holes and some number 
of fc-string with k ^ 2. See Table IHl for a summary of chiral excitations (including also 
the ones discussed below). How many strings of each type can be inserted in a particular 
state depends on in a complicated way and requires a careful analysis of the Bethe 
equations |121IS]- The dependence on N leads to a non-trivial exchange statistics of 
magnons so that the Hilbert space of 2N magnons is not just 2N copies of a single- 
particle space [nUEni- An obvious consistency condition on a distribution of holes and 
fc-strings (but not the only one!) is that the total number of 2-strings, L/2 + M2, should 
be integer and the total number of /-roots should not exceed L + Mu/2. 



6.3 Adding derivatives and spin separation 

Let us finally consider the modes excited by the insertion of a covariant derivative or, in 
the language of the Bethe ansatz, by the addition of a u-root. Composite QCD operators 
associated to physical states of this kind contain one or more derivatives on top of the 
anti-ferromagnetic vacuum of the SU{2)l spin chain. We will see that these states are 
similar to BMN operators in A/" = 4 SYM but with important differences which arise 
because of the anti-ferromagnetic nature of the ground state. The derivative can be 
thought of as an "impurity" propagating on the background composed of /a/j's. Since 
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the anti-ferromagnetic vacuum is a Lorentz scalar, the presence of a derivative adds a 
chiral and an anti-chiral index to the operator. We write: 



^(n) ^ ^ . f D ■ f f w,°o-"2i (6 24) 

^aa J aia2 ■■■ -^aoaj 0121-10121 ■■■ J a2L-iOt2L'r a i V^-'^^) 



1=1 



where we have introduced the wave-function ■?^2o...o2l Qinear combination of products of 
CoiOj) to denote that all but one a^'s are contracted with each other and the remaining 
afc = a. The operator ()6.24|) with zero momentum (n = 0) and ^/^^o-'-^az, = ^cjo^ai---"2z,^ 
where ip^^-'^^L jg ^j-^g ground state wave function, is a descendant: 



^ ] tr /g^Q-g . . . D^fa2i_iOL2i ■ ■ ■ fa2L-ia2L'^gs ~ d^tl fa^a^ . . . /a2L-i"2L V'gs^ (6.25) 

1=1 



and has the same anomalous dimension as the unperturbed operator. In complete anal- 
ogy with the BMN states if the momentum of the derivative insertion is non-zero 
but small, the energy will be small. Of course, the cyclicity of the trace (zero total 
momentum condition) requires taking ip^o-'^^L (different from S^^'ipgs'""^^ ^ i-^- requires 
perturbing the anti-ferromagnetic vacuum. 

Thus, we encounter an important difference between ()6.24|) and BMN operators in 
A/" = 4 SYM. Unlike BMN operators, ()6.24|) contains not one but two elementary excita- 
tions! One of them carries momentum p = 27m/ L and is associated with the dotted index 
carried by the derivative. The other is associated with the extra undotted index induced 
by the derivative, but propagating independently. The reason is that the number of 
linearly independent states obtained by different ways of contracting indices a^, . . . , 
is enormous for large L and in the overwhelming majority of cases the unpaired index 
a that is left after all contractions are done is not the undotted index of the derivative. 
The left spin of the derivative thus looses its individuality and dissolves in the sea of 
the background spin. The excess of the left spin is carried away by an independent 
excitation. The mechanism of spin separation just described is a direct analog of the 
spin-charge separation in [22] • In this particular case the left and right spinors should 
have opposite momenta because of the trace condition, but in a more general setup with 
several derivatives different excitations are absolutely independent. 

We just saw that states consisting of only one derivative must necessarily contain a 
spinon and thus their anomalous dimension scales as 1/L. However, an operator of the 
following type with two derivatives and all undotted indices contracted 

n ■ —\^C, fr f D f D f f J.ai...a2L+2 

aP / , '^Im Jaia2 ■ ■ ■ ^ ol2l+\oi J 021-10121 ■ ■ ■ -^q2I,+2/3'' "^a^-i^Zm • • • Ja2L-ia2L r 

(6.26) 

might have parametrically smaller anomalous dimension ~ since at least the zero 

momentum condition can be satisfied by making the anti-chiral excitations propagate in 
opposite directions. 

To investigate whether it is really possible to have a physical state with energy scaling 
like over the ground state (i.e. without spinous) we study the distortion in the 
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distribution of 2-strings caused by the presence of one w-root exactly in the same spirit 
as in the previous section. 

If we add a w-root, the equations for the density of 2-strings is modified by the 
scattering term. The scattering phase is 

52„(A) = 2 arctan A - tt, (6.27) 

and modifies the equation for the distribution of 2-strings as follows: 

113 1 



2 A2 + 1/4 2 A2 + 9/4 



1 1 

+ 



.(A-O^ + l (A-0' + 4 
The above equation is solved by 



(6.28) 



Pu{X)^p{X) + jCT^{X-u), (6.29) 

where 

is the distortion in the sea of 2-strings caused by the scattering off the u-root. The 
distortion of the density is normalized as 

/+0O I 
dA(7„(A) = -. (6.31) 
-oo ^ 

Thus, the distortion of the Fermi sea caused by the -u-root effectively adds an extra 
1/2 of an /-root to the background distribution and the Lorentz spin of an excitation of 
this type turns out to be (0, 1/2). The scattering off the background distribution of spins 
completely screens the left spin of the derivative, so that the derivative excitation is a 
right spinor, rather than a vector. Note however that this excitation is not viable on its 
own due to the fractional /-root. It needs to be accompanied by a fractional excitation 
of type / or another derivative excitation. A single derivative excitation can only be 
accompanied by a hole in order for the Fermi sea to be occupied by an integer number of 
2-strings. Globally, this state has spin (1/2, 1/2) and dimension L-|- 1 as expected when 
adding a derivative. However, it carries not one, but two independent momenta, one for 
the H-excitation and one for the /-distortion. 



^^In fact the two momenta are related by the momentum constraint, which would otherwise, in the 
case of a single excitation, force the momentum to be zero. 
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We can now compute the energy and the momentum of the derivative excitation: 

Pu[u) = vr — 2 arctan — 

/+00 / 2A\ 

d\au{\ — u) ( TT — arctan 2A — arctan —J 

= 2 arctan h 2 arctan 2u — arctan sinh nu. 

2 3 

/X 3/2 1 /•+°° ,^ X 1 



m2 + 9/4 2 ' V^^ + V4 A2 + 9/4 

+ W—T—, (6.32) 



m2 + 9/4 ■u2_^x/4 2cosh7rM 

which gives a parametric representation of the dispersion relation e = eu{p)- Thus, 
in the hmit of small momenta, the contribution of a right spinor with the momentum 
p = 27m/ L to the anomalous dimension is given by a BMN-like formula: 

- 7o = — . (6.33) 

Looking at the quantum number summarized in Table El we see that the configuration 
with 2 u-modes and a 1-string of /-roots over the anti-ferromagnetic vacuum correspond 
to a physical state with no holes in the 2-string vacuum. Thus there exist operators with 
two derivatives and anomalous dimension of order 0{1/ L?) over the ground state. 

For completeness let us mention that the computation for an fc-string of w-roots is 
similar. We obtain the dispersion relation 

2u 2u 
Pu{u) = 271 — 2 arctan 2(1 — 5k,2) arctan ■ 



k + 2 ' k-2' 
{k + 2)/2 {k-2)/2 

^"^""^ ~ u^ + {k + 2Y/A + u' + ik- 2)2/4 ' ^^'^ ^ 

and find that also one Z-root is effectively added to the Fermi sea. The fc-string describes 
an /c-particle bound state of magnons. 
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A The conformal group and QCD 

The conformal algebra 5u(2, 2) = so(2,4) consists of the Lorentz generators £^ and 

the dilatation operator V, the momentum and the conformal boost /C"^. They obey 
the commutation relations 

[C^,, Ct,] = S'^C', - 6',q, = 0, 4 = 6^cl - s^q, 

[p,/:^] = o, [v,c<^ = o, [v,vj = v,,, [p, /c"^] = -/c"^ 

The conformal algebra has rank three and the Cartan subalgebra can be chosen to 

be 

Jo = {/:}, ^i, V] (A.2) 

(recall that C\ = —C^ whereas raising and lowering operators can be written respectively 

as 

J+ = {Cl Cl /C"^} and J- = {Cl £f, V^^]. (A.3) 



A.l Multiplets 

A highest weight state |h.w.) is a state that is annihilated by the operators in J^. By 
acting with the lowering generators J7L one can generate all the states of the multiplet 
(module) 

{1, J_',...}|h.w.). (A.4) 

We will give an example directly related to QCD below in App. IA.31 

A highest-weight state and thus the corresponding multiplet is characterized by the 
charge eigenvalues of the Cartan subalgebra Jq. A standard way of expressing these are 
the Dynkin labels 

[p,r,g] (A.5) 

where p and q are twice the SU{2)l x SU{2)r spins (5*1, 5*2) and r is a negative number 
related to the scaling dimension D hy r = —D — (p + q')/2. 

A.2 Unitarity and shortening 

Generic unitary representations [p, r, q\ of the conformal group satisfy the unitarity 
bound^^ 

p + r + g<-2, (A. 6) 

^^For the compact counterpart su(4) = so (6), all Dynkin labels would have to be non-negative integers. 
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but there are also exceptional unitary multiplets whose Dynkin labels are related by 

p + r = —l, q = or r + q = —1, p = 0. (A. 7) 

In QCD, the elementary fields correspond to the second type while local operators are of 
the first type. Note that all the unitary multiplets are infinite-dimensional as required 
by the non-compact nature of the conformal group. 

Generically, all combinations of the momentum generator V generate new states when 
acting on the highest-weight state. Here, there are four exceptions 



scalars 


P 











r = 


-1, 




chiral fields 


P 











r = 


-1 - 


P, 


anti-chiral fields 


P 











r = 


-1 - 




conserved currents 


P 











r = 


-2 - 


q-p. 



We know for instance that the scalars^'' satisfy the equations of motion (classically) 
P^lh.w.) = while the conserved currents are defined by P'^jh.w.^) = 0. The chiral and 
anti-chiral fields satisfy both types of conditions. 

The chiral and anti-chiral fields ()A.8jl are the central objects in QCD. We will denote 
these multiplets by 

: [A;,-l - A;,0] and V"^ : [0, -1 - fc, A;], (A.9) 

where k is restricted to 1 or 2 in the specific case of QCD. The fields also enjoy a number 
of useful features in terms of representation theory which makes them rather easy to 
handle despite the fact that they are infinite-dimensional. For instance, they can be 
represented by means of a set of harmonic oscillators as we shall demonstrate below. 

A. 3 Oscillator representation 

The conformal algebra has a nice representation in terms of a set of bosonic oscillators 
[a", ajj] = 5| and [6", 6p ~ = 1, 2, d = i, 2). The generators take the form: 

D = 1 + lay + ifetfei, P^^ = at6t, /C"'^ = a V (A.IO) 

which can easily be seen to obey the commutation relations ()A.1|) . (See e.g. for the 
general theory and extension to supergroups.) It is also possible to construct an extra 
u(l) generator, not part of the conformal algebra and commuting with it: 

A = la{a^ - \b\b\ (A.ll) 

This generator is associated to the "chiral" type symmetry transformation discussed in 
the text. 

^^Of course, there are no elementary scalar fields in QCD but it is just as easy to discuss the general 
case. 
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If we choose a vacuum vector |0) annihilated by a"" and b^, highest weights are of 
type (c!'2)'^|0) However, since the Lorentz structure is the famihar one, we 

will only use /C°^|h.w.) = as the relevant condition and refer to the whole Lorentz 
multiplets aj^^ . . -(lajO) and ftj^^ . . -bl^jO) as "primaries". The modules constructed by 
acting with lowering operators on the above states provide the oscillator realization of 
the fundamental fields in the Lagrangian. Thus, the modules ()A.9|) are spanned by: 

= {al . . . aijO), al... a^J^^ |0), . . . al^/,/,^ |0), . . . >, 
^-' = {bl...bl\0), albl...bljO), alalbl---bljO),...). (A.12) 

For convenience, the modules appearing in QCD shall also be denoted as V* where 
^ = fi I iXii^iXi^ are the physical primary fields. 

Recalling that Pa = —iD^n we write (e.g. fap = o-aOfllO)): 





= = 




y/ 


= = 


(^faia2i ^ f 01611612013: D f a.\a2Ci\ci26iz6n 




= = 


{V^Qi ) Dlp(^-^(^2ai ; D 'ipaia2a3aia2 ' ' ' } 




= = 


{Xai 7 ^X,aia2ai j D Xaia2a36ii6i2 ■ ■ ■ } 




= = 






= = 


{V^ai5 -0^«i«i«2 5 ^ ^0102010203 ■ ■ ■ } 



The indices of the same kind (dotted or undotted) in the above modules are always 
understood as totally symmetrized. 



A. 4 Tensor products of two fields 

In computing the action of the dilatation operator on two adjacent fields in a gauge 
invariant operator we are naturally led to the problem of decomposing the product of 
two of the modules ()A.13|) . In the oscillator representation we would have to introduce 
two commuting sets of oscillators [c^fi)) ctjj) /?] = ^ij^'/s and ^] = ^ij^"^ where we 

will only be interested in the case i,j = 1,2 but in principle one could generalize to 
the product of an arbitrary number of ()A.13|1 . The generators are simply the sums of 
the generators for each set, in particular the conformal boost for the product of two 
singletons ()A.13|) is JC"^ = + (^'(2)b(2)- Highest weight states in the product can 

be found by looking for solutions to /C"^|h.w.) = amongst all linear combinations of 
states of the type 

«(l),o, • • • «(l),o,&(l),oi • • • ^(l),o,«(2) A • • • «(2),/3™^[2) A • • • ^(2),/^JO) ^^'^"^^ 

where now |0) is annihilated by both sets of a^^^ and 6"-^. 

Obviously there are more possibilities to satisfy the highest weight condition and in 
fact the product of two arbitrary modules ()A.13|) decomposes into an infinite sum of 
irreducible components, distinguished by a quantum number j named "conformal spin" . 
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After using the oscillator representation to get some practice, one can recast the results 
in the more compact language of fields or in the even more compact language of Dynkin 
labels. Assuming k > \k'\ we have 

CO 

V±^® V'^' = ^ V^^^'+^''. (A.15) 

j=Tk' 

The Dynkin labels of the modules in the decomposition ()A.15|1 are 



[s + 2j,-s-2-j,0] for-s/2<j<0, 
[s+j,-s-2-2j,j] forj>0, 



, ( [0,-s-2-j,s + 2j1 for-s/2<j<0, , 
' [[j,-s~2-2j,s+j] forj>0. 

A proof of these identities is sketched in App. IA.51 

Let us explain in practice how the decomposition works by considering the product 
of two modules V-^: 

oo 

i=-2 

To understand the structure of the decomposition we present the explicit form of the 
first few elements in the first few modules. In order to do that we shall introduce the 
following notation: 

• All indices of the same kind written explicitly are meant as totally symmetrized 
even if they belong to different fields. 

• Anti-symmetrization between two indices by contraction with e"^ is denoted by 
replacing the contacted indices by the symbol x (or x). 

For instance e^^ f.y(af(3)5 is written as fxaf/sx- As long as we consider only the products 
of two modules ()A.13jl . there is no room for confusion since indices can only be anti- 
symmetrized between distinct fields. 
Thus we have: 

fxx+ fxxDfxxaa, • • • ) , 

Dfxxaafxx " fxxDfxxaa, ■ ■ ■ ) , 
~ {falsf^S, D faf^nrif-fS + f ajsD f^Srir)^ D fxajSrifx-y ~ f xcxD f x p-^r)-, ■■■)■, 
={Dfal3r}r)fi5 — fa/sD f^Srifj, ■ ■ ■) (A. 18) 

and so on. The first element in each module is the primary. The descendants are obtained 
taking derivatives and decomposing the result into irreducible representations. Primaries 
in the modules Y^^ for j > have j covariant derivatives. Similar decompositions hold 
for all the modules. In particular: 



oo 



j=+2 
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where, for instance, 

W = {Uhp. DU,J^. + UDf^^^., ...) (A.20) 

and so on. 

Similarly, we shall write V-^ = Yl'jL-i^j^ so on for the tensor products 
involving fermions. The lower bound of the sum can be read off from ()A.13IA.15|) in 
each case. 



A. 5 Polya counting 

Using Polya counting^^ we can prove the decomposition of tensor products ()A.15j) . Let 
us denote a state with weight [p, r, q] of SU{2, 2) by the monomial a^b'^d^^^^^'^^''^'^ 

[p, r, q] oPWd-'"-^'^-''!'^. (A.21) 

In other words the exponents of a, h are the third components of spin and the exponent of 
d is the dimension. A multiplet of states is consequently written as a polynomial, where 
the integer coefficients indicate the multiplicities of states with fixed quantum numbers. 
This notation is very convenient for dealing with the infinite dimensional modules that 
appear for the conformal group. For example, consider the field strength component /n 
and the set of all descendants using the derivative 

{/n,Dn/n,/^n/^n/ii,...}. (A.22) 
Using the rule ()A.21jl we can write this as the polynomial 

a^d^ + aHd^ + a%^d'' + . . . . (A.23) 
Of course, this is just a geometric series which sums up as 

c?d^ 



dab 



(A.24) 



So we have found a compact way to represent /n along with its D^i descendants. 

All the states of a long unitary multiplet [p, r, q] of the conformal group are summa- 
rized in 

^'""'^^ ~ {1 - dab){l - da/b){l - db/a){l - d/ab){l - a-^){l - b-^) ' ^ '^^^ 

In the denominator one finds the momentum and SU{2)l and SU{2)ii ladder generators. 
The two differences in the numerator make SU{2) multiplets finite-dimensional. 

Generic unitary multiplets [p, r, q] of the conformal group satisfy the unitarity bound 
flA.6|) but there are also the exceptional multiplets in ()A.8|) which need special attention. 
Let us discuss these special cases, where we refer to the long multiplet defined in ()A.25j) 

""^^See e.g. IHSI for an introduction in the context of field theory. Polya theory has also been used in ]^ 
to count single-trace operators in A/" = 4 SYM. 
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as [p, T, g]L. The simplest example (not directly related to QCD) is an on-sliell scalar, 
which has Dynkin labels [0, —1, 0]. It turns out that in order to remove the terms that are 
reducible due to the equation of motion D'^cp = ... we have to subtract a long multiplet 
[0, -3, 0]l.^^ 

[0,-1,0] = [0,-1, 0]l- [0,-3, 0]l. (A.26) 
For an on-shell chiral field of spin 5*1 = p/2 with r = ~1 — p one finds 

[p, r, 0] = [p, r, 0]l - [p - 1, r - 1, 1]l + [p - 2, r - 1, 0]^ (A.27) 

and analogously for an anti-chiral field. These multiplets violate the bound ()A.6|1 . but 
due to their tightly restricted charges they are indeed unitary. 

A less restricted type of multiplet, a conserved current, sits right at the unitarity 
bound ()A.6|) . Its state content is reduced by the states D'^J^... = 

[p, r, q] = [p, r, q]^ - [p - 1, r, g - 1]l. (A.28) 

Using these expressions and sum rules of geometric series, it is straightforward to 
prove ()A.15|1 . 



B Length-two primary states 

Let us now discuss the gauge invariant primary states that can be constructed with only 
two fundamental fields (L = 2) and an arbitrary number of derivatives. We do this 
partly to present a check of our results against the standard literature and partly as an 
illustration of how to use the Hamiltonians ()2.12j) and ()2.13p . 

Reducing to L = 2 chains amounts to projecting all modules in the decompositions 
()A.16|) on the gauge singlet states. In considering L = 2 states and in order to make 
connection with the classical papers on the subject is more convenient to rewrite the 
gluon modules in a way that makes parity manifest. To this end, we introduce the 
following projection operators: 

p//±/7 = i (I ± (_)^x) (Pf + Pf ), (B.l) 

in terms of which the Hamiltonian ()4.4|) reads 

oo oo 

H^^ = {P" + Pf ) + E [eM^'^" + E-Pf-") , (B.2) 

j=-2 j=+2 



where 
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Not incidentally this is just the weight of 
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For sake of clarity, we give expression for the primary and the first few descendants of 
these new modules: 



V: 



■//+// 



V: 



//-// 



y//+// 



{faplap + /a/3 /a/3 5 

-^/a/37x/xd ~ fxa^fajS-yx^ faxD f^^0^ — D f^^0^fax ^ • • • ); 

DfaP^afp^ + fapDf^^p^ — /a/j-D /a/377 " ^faafi^fpii 
^faf^'yxf-ka + fxa^faf^'yx^ fax^fxa$^ ~^ ^fxa$^faxj • • • ); 

{Dfaf3-rafpy ~ foipDj^^p^ + /q,/3-D /a/377 ~ /aa/S-y //37 ' • • • ) ; 



where the conventions about symmetrization and anti-symmetrization are the same as 
those explained in Appendix 1X1 

For the gluon modules the projection onto gauge singlet states amounts to taking the 
trace. All the elements in are traceless (thus projected out) for j even and those 

in Vj"'^, Yj^ and Y^j^^^^ are traceless for j odd. Of course, the quark-gluon modules 
contain no gauge singlet at all and the quark-quark modules contain singlets for all j's. 
The surviving length- two gauge invariant primaries are summarized in Table [7| 



Name 


Range 


(^1,^2) 


D 


[P, r, q] 


r 


X 


A 








(0,0) 


4 


[0,-4,0] 


4 





2 


¥/ 


3 


= 0,2,4... 


(2+j72,j72) 


4 + j 


[4 + j, -6-2j,j] 


2 


4 


2 








(0,0) 


4 


[0,-4,0] 


4 





-2 




3 


= 0,2,4... 


(j72,2 + j72) 


4 + J 


[j,-6-2j,4 + j] 


2 


-4 


-2 


yff+ff 
j . _ 
y//-// 


3 


= 2,4,6... 


(3/2,3/2) 


2 + 3 


[j, -2-2j,j] 


2 








3 


= 3,5,7... 


(3/2,3/2) 


2 + 3 


[j, -2-2j,j] 


2 














(0,0) 


3 


[0,-3,0] 


3 





1 


yX^ 


3 


= 0,1,2... 


(1+j72,j72) 


3 + J 


[2+j, -4-2j,j] 


2 


2 


1 






(0,0) 


3 


[0,-3,0] 


3 





-1 


i 


3 


= 0,1,2... 


(j72,1+j72) 


3 + J 


[j, -4-2j,2+j] 


2 


-2 


-1 


i _ 


3 


= 1,2,3... 


(3/2,3/2) 


2 + 3 


[j, -2-2j,j] 


2 








yXX 


3 


= 1,2,3... 


(3/2,3/2) 


2 + 3 


[j, -2-2j,j] 


2 









Table 7: Table of all gauge singlets primaries of length two. The indices (^i, S2) refer to the 
Lorentz spins, D is the classical dimension, [p,r,q] are the Dynkin labels and A is the U{1) 
charge introduced above. The twist is defined, as usual, hy t = D — {p + q)/2. In the presence 
of "asymmetric" operators, dimension and twist are not enough to uniquely specify the Lorentz 
representation - we need an extra quantum number, e.g. the "asymmetry" factor X = P ~ Q- 
Notice the presence of exceptional cases with negative conformal spin. 
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We are particularly interested in computing the one-loop anomalous dimension of the 
primaries. First of all, in order to disentangle the contribution of descendants carrying 
the same quantum numbers, one always performs the computation at zero injected mo- 
mentum, for which total derivative terms simply do not contribute. Because of the U{1) 
symmetry relating the anomalous dimensions of the related operators, we can restrict 
our attention to A > 0. Finally, the anomalous dimensions of the primaries in the last 
two modules in Table [7| are obviously the same. We are thus led to consider the renor- 
malization of the following operators at zero injected momentum:^'' (the range of the 
conformal spin j is the same as in Table [7j) 

V{{ : tr F^^F^, 

Vf : tr F^^.'Op, . . . ^ p^F^)„ - traces 

Vf+^"^ : tr -traces 

: tr -traces 

V^t : (1(1 

VJ'^ : g[7A., l{v] D p^... D p.)q - traces 

Vf : g7{p,'^p,..-'^p,)g- traces. (B.4) 

Indices inside round brackets are totally symmetrized and by "—traces" we mean that the 
tensor structure is irreducible, i.e. tracing over any two indices yields zero. This reduction 
will always be understood in the following. The operator tr F^^F^y (the Lagrangian) 
contains both the selfdual and anti-selfdual part but both terms renormalize the same 
way due to the f/(l) symmetry. Its anomalous dimension is well known: 

tr F^^F.^ : 7 = ^ = 2^^E^{. B.5 

This is the first check. The factor of two comes because we must sum the (equal) 
contributions of Hf| and Hf^. Similarly, the quark mass term qq contains both the 
chiral (xV^) and anti-chiral {x4') and its anomalous dimension is also well known: 

(1(1- 7 = -9^7— = (B-6) 
Z Zn Ztt 

(This time there is no factor of two because we are considering an open chain.) 

The anomalous dimensions of the primary operators of Yj^~^^^ and V^^ are also well 
known from the classical work of pp. We do not write explicitly the flavor structure 
but even for flavor singlets, in the large Nc limit, the mixing between gluonic and quark 
operators is suppressed. In particular, in this limit their matrix of anomalous dimension 
is upper triangular and the eigenvalues can be read off from the diagonal entries. It is 
even possible to rescale the quark kinetic energy in the Lagrangian by a factor \/7^ and 

^"^The exact conformal operators contain total derivatives in particular combinations fixed by confor- 
mal symmetry |37|. 
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obtain a diagonal matrix: 



^7(pi D p^... D p^)q : 

7 = ^ (Mj + 1) + Mj - 1) - 3/2) = ^ £;f . (B.7) 
Ztx ztt 

In the same way, the anomalous dimension of the primary operator of can be 

read off from the work of |57j : 

tr Fl^^Vp,...Vp^_,Fp^)p: 



7 = ^^ r + 4/i7 ^2^^^-. B.8 



The work [57, contains also the matrix of anomalous dimensions of a second group of 
operators but the gluonic term appears to be a descendant, casting some doubts on the 
relevance of this second set of anomalous dimensions. 

We are left with the anomalous dimensions of the primaries for and . We 
worked out their anomalous dimensions by diagrammatic techniques, employing the pre- 
vious trick of zero injected momentum and the following additional ones (also well known, 
see e.g. [58 ): 

• The symmetrization of the indices was accomplished by contracting each of them 
with a constant vector A„, etc. 

• The tracelessness condition was enforced by dropping in the evaluation of the 
Feynman diagrams all terms containing (A)^, (and for the second case also 

• Finally, to ensure that the anti-symmetrization of any three indices had been re- 
moved, we evaluated a particular component that is manifestly free of such con- 
tribution, namely we set (at the end!): Ap = 5ip, fj, = 2, (and for the second case 
also (T = 2). 

The resulting anomalous dimensions are: 

7(.]K • • • Dp,g : 7 = ^ {^hU + 1) - ^ ^ Ef, 
F^^Jd p, . . ."d p^F,)^ : 7 = ^ + 2) - ii") = 2^ , (B.9) 



27r V 3 y 27r 

confirming the corresponding entries in the Hamiltonian. The last two operators are 
"asymmetric". They cannot be generated perturbatively in the operator product of 
two currents because such product does not change chirality. They might however be 
interesting for comparing with lattice calculations 
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C R- matrix 



Formalism. Let us give a brief introduction into the R-matrix formalism. We define 
several vector spaces which we label by k,l, . . .. These are (not necessarily equal) modules 
yfc,;,... q£ ^ symmetry algebra and they are thus called 'spins'. In our case the symmetry 
algebra is the conformal algebra and a 'spin' will be a chiral field strength multiplet. To 
each spin we associate a spectral parameter Uk^i^...- The scattering of two spins k,l is 
assumed to be elastic: Neither the spin modules V*^'' nor the spectral parameters Uk^i are 
modified, there is merely a phase-shift when two spins are interchanged. The phase shift 
is described by the R-matrix "R^kiiuk — ui), which is a unitary bi-linear operator acting 
on the two spins 

B.ki{uk - ui) -.Y'xY^Y'x Y. (C.l) 

In this formalism, integrability means that the Yang-Baxter equation is satisfied. It 
ensures that the order in which spins scatter does not matter 

where we have defined 

Uki = Uk- ui. (C.3) 

The form of the R-matrix is in general very difficult, but it can be constructed for 
arbitrary semi-simple algebras and arbitrary irreducible representations. In our case the 
'spins' will all be of the type of oscillator representations. These have special properties 
that allow us to compute the R-matrix rather conveniently. 

R-matrix for spin oscillators. First of all, we introduce the compact notation JT^b 
for the generators of su(2, 2). Here, an index A combines a chiral and an anti-chiral index 
a and a. We can now introduce a combined oscillator A-^ = {a", 6^^} which behaves like 
a four- component version of a°. The combined generators are now written as 

J^B = A^B^^ - \ A = \A\A^ + 1 (C.4) 

and satisfy the commutation relation 

[J^B^ J'^d] = S^J^'b - S'iJ^D, [J^B, A] = 0. (C.5) 

Let us now consider the R-matrix acting on the two oscillator modules and WK We 
know that the R-matrix is invariant under conformal symmetry therefore it is useful to 
consider the decomposition of the tensor product V'^ (S) V into irreducible components, 
flA.15j) . As each component V^' appears with multiplicity one, the R-matrix can merely 
assign one eigenvalue to each: 

^kiiu) = EjPfRfi^)- (C.6) 

We will show below that these eigenvalues obey the recursion formula 

{Ukl) - _i(j2 _ j2 ^Rj^l[Ukl)■ ((-.V) 

"-kl 4 ['JkLj •^kLj-l) 
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Here J^i j denotes the eigenvalue of the quadratic Casimir 

J' = J'^bJ^'a (C.8) 
on Vj'. For a highest-weight representation with Dynkin labels [p,r,q] it is given by^^ 
4,r,,] = Hp + 2) + i (2r + p + g) (2r + p + g + 8) + + 2) . (C.9) 

This result for the R-matrix is completely general, it holds for any two oscillator repre- 
sentations, even fermionic oscillators, for any unitary group. 

Closed Chiral Chain. Let us now consider only spins which are chiral field strengths, 
i.e. only modules of the type V-'^ = V[2,-3,o]- The monodromy matrix of this system is 
given by 

n^iu) = Rfc,i(u)Rfc,2(M) • • . Rk,M (C.IO) 

where the auxiliary spin labeled by k transforms like a field strength . The transfer 
matrix is 

T(m) =trfcfifc(M) (C.ll) 

The nearest-neighbor Hamiltonian of an integrable system is given by the logarithmic 
derivative of the R-matrix at m = 

H// = zR7; (0)^(0). (C.12) 
The recursion formula for the eigenvalues Ej^ of the two-site Hamiltonian is thus 

= ,2 ^ p + E^U- (C.13) 

In our case the quadratic Casimir ()(I9jl of V|-^ is given by J| = 2(j + 2)(j + 3). The 
recursion formula for the energy eigenvalues is thus 

= -^ + E"=2h{j + 2) + 2Ef. (C.14) 
J + 2 

where E-l" is an energy shift that does not follow from integrability. This agrees precisely 
with the result from planar QCD and thus the sector of chiral field strengths is integrable. 

Open Chiral Chain. The integrable open chiral chain has spins Y^ at all sites except 
1 and L, where the spin is = Y^ = V[i^_2,o]- This system is characterized by the 
monodromy matrix f^^^ 



n^iu) = Rfc,i(M ± |)Rfc,2(M) . . . Rk,L^^iu)RkAu ± i) (C.15) 



^"'^A curious observation is that the quadratic Casimir vanishes precisely for [0,0,0] (trivial represen- 
tation), [2,-3,0], [0,-3,2] (field strength components), [1,-4,1] (conserved flavor currents), [0,-4,0] 
(Lagrangian) and [1, —2, 1]. 

^^In 0H| the system of SU{2) spin +1 in the bulk and SU{2) spin +1/2 at the boundary was investi- 
gated. It arises in the SU{2) l sector of QCD and the SU{2, 2) system is a straightforward generalization. 

^■^For the chain of length L = 2 we should choose ^k{u) = Rfe,i(u)Rfc,2(u) instead where the auxihary 
spin is = V'''. 
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where the auxihary spin labeled by k transforms like a field strength V-^ . The open chain 
transfer matrix with SU (2, 2)-preserving boundary condition (K = 1) is 



Tiu) = tiknkiu)n^\-u). (C.16) 

The Hamiltonian in the bulk is given by the same expression as before. At the boundaries 
the Hamiltonian is modified due to the shift in the definition of Qk{u) 

H,^ = ^R,J{±^/2r'^{±^/2) (C.17) 
consequently, the energy eigenvalues obey 

The eigenvalue of the second Casimir for is J'^j j = + 8j + when we substitute 
this we obtain 

Ef = + + Efl^ = h{j + 2) + hU + 1) + E^ + Ef (C.19) 

J + 1 J + 2 ■' 

where E^ + E-^ is a constant we cannot determine here. 

Also the chains with one or two anti-chiral fermions at the ends, but all chiral field 
strengths are integrable. We simply replace first/last R-matrix in Q by Rfc,i/L(M ± f-)- 
The eigenvalue of the second Casimir for Yf is r • = 2j^ + 4j — | and we obtain 



E: 



- + + Efl^ = h{j + 2) + h{j -1) + E^ + Ef. (C.20) 

J + 2 J - 1 



In this case an anti-chiral field is not an obstacle to integrability. The reason why an 
anti-chiral field strength breaks integrability is that it may propagate while the anti- 
chiral quark is fixed at the end of the chain. In other words, integrability breaks down 
for mixed field strengths due to E^-^'^ ^ 0. 

Proof. Here we will present a proof of the Yang-Baxter equation for the R-matrix 
given above. It is similar to the one for J\f = 4 given in JH]- The three spin modules 
for the YBE will be of oscillator type and one of them will be fundamental. First of 
all, consider the R-matrix of an oscillator spin k and a fundamental spin i. The tensor 
product V'^ ® Vfund has two irreducible components which we shall denote by V^*. The 
projectors are given by 

2 2y4/j — 1 

where Ak is the value of A on V'^ and J'ki is the algebra generator acting on whose 
components jTfc^s specify the map Vfund — ^ Vfund- The difference of eigenvalues of the 
quadratic Casimir is 

4,+ - Jl = 2{2A, - 1). (C.22) 
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The R-matrix ()(I6|(I7jl therefore reads 

Uki + ^[.J ki,+ ~ J ki,-) Uki + i^y^Ak - I) 

Now consider the YBE for two oscillator spins fc, / and a fundamental spin i and 
expand the R-matrices Rjfc(MiA:), R-iK'^iz) involving a fundamental spin 

= 'R,ik{uik)^ii{uii)'R.ki{uki) — Hki{uki)^ii{uii)'Rik{uik) 
_ Ouik + Uii + 2i- ^Ak - ^Ai) 



2{u,k + U2Ak-l)){uu + ^{2Ai-l)) 



\Jki + Jlh '^kliUkl)\ 



2{u,k + ^{2Ak 


-l)){ua + ^{2Ai 
1 


-1)) 


2{uik + ^{2Ak 


~l)){uu+'^{2Ai 
1 


-1)) 


2{uik + {{2Ak 


-l)){ua + {{2Ai 


-1)) 



{[iTfcj, iTii], 'Kkiiuki)^ 
[{Jki.Ju},'£lki{uki)\. (C.24) 
The first term vanishes because 

Jki + Jli = Jkl,i (C.25) 

is the rotation generator acting on ® which clearly commutes with the covariant 
R-matrix. Let us now introduce the symbol Cki,i for the combination of generators in 
the second term 

Jki — Jli = Cki^i- (C.26) 
Then we can express the commutator in the third term in terms of Cki^i 

[Jki, Jli] = zlJkhCki^i] (C.27) 

and the quadratic Casimir J^i on xYK Finally, we can write the anti-commutator in 
the fourth term using squared generators'^ 

{Jki, Jli} = {Jkl,iY — {Jk,iT ~ iJi,i)^- (C.28) 

Using the rules ()C.4IC.5|) for oscillator representations, we can compute the contracted 
product of two generators 

J^bJ^'c = iA- 2)J^c + IM' - l)^c- (C.29) 

It is a special property of oscillator representations that J^bJ^c can be written as a 
linear combination of J'^c and 5^.'^ We can now simplify the anti-commutator ()C.28|) 

{Jki, Jli] = {Jki^i? - {Ak - 2) Jfc,, - (A - 2)Ji^, - 1(^2 + ^2 - 2) (C.30) 
= {Jki,if - \{Ak + Ai- A)Jki,i - |(v42 + Af-2)- liAk - Ai)Cki,i. 



^"^Note that we have to distinguish between (Ji)^ = JiJi and the quadratic Casimir = tr^ JiJi. 
25pQj- generic representations only a four-fold product of generators of the above sort can be written 
as something simpler. 
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Again all terms except the one proportional to Cki,i drop out in commutators with Hm- 
We can now put everything together and find 



— Tiif:{uik)tlii{uii)'Rkc{uki) — 'Rki{uki)Ilu{uii)'Rik{uik) 

_ . Ukl [Ckl,i, HkliUkl)] + I [j'kh {Ckl,i, ^kliUkl)}] 

2{uki + |(2A, - 1)) {uu + f (2A - 1)) 



(C.31) 



To proceed further, let us investigate the form of the operator Cki^i- It is a map (S> 
V' ® V* — s> ® V V*. By construction it is invariant under simultaneous rotation 
of all modules by the same amount. When we combine the (oscillator) modules into 
irreducible components, we can write 

Cki, : Vf ® Vf,„d ^ Ef^r ® ^fnnd- (C.32) 

By comparing the highest weight vectors of V^' and V^/ we see that Cki,i can only be 
invariant when j' G {j — + 1}. Furthermore, Cki,i has negative parity under 

interchange of k and /. As the parity is alternating with j, the cases j' = j ± 1 are 
allowed while j' = j is not. In conclusion we can write 

Ck^ : Yf ® Vf,nd ^ Yf_, ® © Vj^i ® Vf,nd (C.33) 

This is another generic feature of tensor products of oscillator representations. The 
recursion relation ()C.7|1 

pkl(^, N _ + \ {Jkl,j - Jkl,j-l} rjkl N nA\ 

Rj [Ukl) - -JY2 f-2 — Y Rj^Auki) (C.34) 

follows from projecting ()C.31|) to the various irreducible components. Similarly, for a 
different ordering of the three spins we obtain the same relation from the YBE 

= Rfcj(Mfcj)Rfc/(tXfc/)R-ii('"jO ~ '^iiiuii)'R,ki{uki)'^kiiuki) 

i{uii — Ukl — ^Ai + ^Ak) r 1 

^ 2(«.. + t(2A,-l))(«, + t(2A-l)) + ^''^""''^^ 
I z{uki + 2z-lAk-lAi) _-T R r./ ^1 

+ 



2{uk^ + ^{2Ak 


-l)){uu + ^{2Ai- 
1 


1)) 


2{uki + ^{2Ak 


-l)){uu + ^i2Ai- 
1 


D) 


2{uk^ + ^{2Ak 


-l)){uu + ^i2Ai- 
1 


1)) 



2(n,., + f (2Afc - 1)) {uu + W^Ai -1))^' 
Uki[Cki,i,'Rki{uki)] + ^{[J'ki, Cki,i],Ilki{uki)} 
2{uki + \{2Ak - 1)) [u,i + \(2Ai - 1)) 



(C.35) 
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